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-  accurate  calibrations  of  peak  prosiQu.ro  (P)  vs  flap  length  (x) 

are  presented  for  the  NOL  large  scale  gap  teat  for  both  pontoi  1  to  and 
tetryx  donor  chargee.  These  calibx'ations  are  derived  from  particle  voloc 
ity  and  ©hock  velocity  measurements  in  the  polymothyl  methracrylato 
(PMMA)  attenuator.  Particle  velocities  are  recorded  by  tho  eloctromag- 
netic  velocity  gage;  shock  velocitiea  are  determined  from  streak  camera 
records  of  tho  progress  of  shocks  through  the  attenuators,  Tho  new 
calibrations  will  allow  direct  comparison  on  a  peak  pressure  basis  of  the 
results  of  teats  performed  in  the  past,  mostly  with  tetryl  donor  charges, 
and  those  which  have  been  performed  recently  with  pentolite  donor  charges 
All  future  gap  tests  will  use  pentolite  donor  charges. 

As  indicated  above,  wo  measured  both  shock  velocity  (U)  and  particle 
velocity  (u).  However,  tho  Hugoniot  of  PMMA  was  selected  by  comparing 
these  U-u  data  with  the  low  pressure  U«u  data  recently  published  by 
Sandla  Corporation;  the  experimental  design  and  accuracy  of  the  latto.r 
work  is  superior  to  those  of  the  gap  test  for  Hugoniot  determinations. 

The  Hugoniot  so  obtained  was  used  in  conjunction  with  our  u-x  data  to 
produce  the  desired  calibration  curves  (P-x). 

The  new  tetryl  calibration  has  about  ?.0%  greater  P  at  140  cards  than 
that  Of  the  previously  published  calibration.  This  new  curve  does  not 
cross  the  old  curve;  therefore  the  relative  ranking  of  explosives  by  gap 
sensitivity  will  not  be  changed  by  tho  use  of  the  new  rotation. 
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CALIBRATION  OF  THE  NOL  LARGE  SCALE  GAP  TEST;  HUGONIOT  DATA  FOR 

POLYMETHYL  METHACRYLATE 


Prepared  byj 

o,  o,  Erkrnan,  D.  J.  Edwards,  A.  R.  Clairmont,  Jr„  ,  and  Donna  Price 


ABSTRACT!  More  accurate  calibrations  of  peak  pressure  (P)  vs  gap 
length  (x)  are  presented  for  the  NOL  large  scale  gap  test  for  both 
pentolite  and  totryl  donor  charges.  These  calibrations  are  derived 
from  particle  velocity  and  shock  velocity  measurements  in  the  poly¬ 
methyl  mothracrylato  (PMMA)  attenuator.  Particle  velocities  are 
recorded  by  the  electromagnetic  velocity  gage;  shook  velocities  are 
determined  from  streak  camera  records  of  the  progress  of  shocks 
through  the  attenuators.  The  new  calibrations  will  allow  direct 
comparison  on  a  peak  pressure  basis  of  the  results  of  tests  performed 
in  the  past,  mostly  with  tetryl  donor  charges,  and  those*  which  have 
been  performed  recently  with  pentolite  donor  charges.  All  future  gap 
tests  will  use  pentolite  donor  charges. 

As  indicated  above,  we  measured  both  shock  velocity  (U)  and  particle 
velocity  (u)  .  However,  the  Hugoniot  of  Pf-5MA  was  selected  by  comparing 
these  U-u  data  with  low  pressure  U-u  data  recently  published  by 
Sandia  Corporation;  the  experimental  design  and  accuracy  of  the  latter 
work  is  superior  to  those  of  the  gap  test  for  Hugoniot  determination. 
The  Hugoniot  so  obtained  was  used  in  conjunction  with  our  u-x  data 
to  produce  the  desired  calibration  curves  (P-x). 

The  new  tetryl  calibration  has  about  20%  greater  P  at  140  cards  than 
that  of  the  previously  published  calibration.  This  is  largely  due 
to  the  much  improved  definition  of  the  wave  interaction  in  the 
attenuator  at  a  distance  of  about  35  mm  from  the  donor  charge  which 
results  in  an  abrupt  change  in  the  slope  of  the  u-x  curve  at  that 
distance.  This  phenomenon  is  also  observed  in  the  streak  camera 
records  of  the  progress  of  shocks  through  the  attenuator;  the  shock 
velocity  vs  distance  curve  derived  from  the  records  has  a  hump  at 
30-40  mm  from  the  donor.  The  new  calibration  curve  does  not  cross 
the  old  curve;  therefore  the  relative  ranking  of  explosives  by  gap 
sensitivity  will  not  be  changed  by  the  use  of  tho  new  relation. 
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The  pen to 1 it©  calibration  curve  gives  larger  values  of  P  than  does 

the  tetryl  calibration  for  x  <  200  cards  (50  mm) .  The  difference  is 

not  significant  experimentally  because  it  is  the  order  of  magnitude 

of  errors  in  both  calibrations.  Excluding  the  nominal  range 

(x  <  10  mm),  the  greatest  difference  is  2.6  kbar,  6r  7%,  at  140- 

150  cards.  For  x  >  200  cards,  the  two  new  calibrations  are  coincident.  » 
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The  work  described  in  this  report  was  carried  out  under  Task  IR-59, 
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Research  Program.  The  results  are  more  accurate  calibrations  of  the 
NOL  large  scale  gap  test  for  both  tetryl  and  pentolite  donor  charges 
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1.  INTRODUCTION 


Two  major  and  related  objectives  of  this  work  are  (a)  to  g  >  ? 4 

recalibrate  the  large  scale  gap  test  (LSGT)  and  (b)  to  select  the 
best  Hugoniot  for  polymethyl  methacrvlate  (PMMA).  A  calibration  ,  ';>■/■ ie 
(or  a  check  of  it)  has  to  be  made  whenever  a  new  lot  of  donor  pellets 
is  received.  The  present  recalibration  differs  from  such  a  routine 
check  in  that  we  are  able  for  the  first  time  to  measure  particle  : 
velocity  directly  and  accurately  (to  ±  0.03  mm/  sec),  and  also  because 
exact  data  (shock  velocity,  U,  vs  particle  velocity,  u)  defining  the 
Hugoniot  of  PMMA  in  the  lower  pressure  range  have  been  published 
only  recently.  Use  of  the  new  method,  the  new  data,  and  improved 
optics  and  data  reduction  for  determinations  of  U  showed  that  both 
of  our  old  calibrations!  (tetryl  and  pentolite  loaded  PMMA)  are  less 
accurate  than  is  justified  by  our  present  improved  capabilities. 


An  accurate  Hugoniot  (shock  velocity,  U,  vs  particle  velocity,  u) 
for  PMMA,  the  attenuator  in  the  LSGT,  is  necessary  to  carry  out  a 
calibration  by  measuring  only  one  of  these  variables  as  a  function  of 
attenuator  thickness.  Unfortunately,  the  Hugoniots  proposed  for  PMMA 
are  legion.  Heretofor  we  have  used  our  U-u  values1  derived  from 
indirect  measurements;  those  values  were  definitelv  less  precise  than 
our  present  ones.  Publication  of  exact  data  in  the  lower  pressure 
region,  in  combination  with  our  own  improved  measurements  of  U  and  u, 
allows  us  to  construct  what  we  consider  an  accurate  PMMA  Hugoniot 
over  the  pressure  range  of  the  gap  test. 


2.  METHODS  FOR  CALIBRATING  GAP  TESTS 

This  section  gives  a  brief  review  of  methods  that  have  been  used, 
or  could  be  used,  to  calibrate  gap  tests.  For  some  of  these,  an 
analysis  is  given  of  the  error  inherent  in  the  method.  One  method 
is  direct:  manganin  pressure  gages  could  be  used. 2  Bv  doing  exper¬ 
iments  in  which  the  gage  was  at  different  distances  from  the  explo¬ 
sive,  a  P  vs  x  relation  would  result.  This  direct  method  is  seldom 
used,  possibly  because  it  is  still  relatively  new.  For  this  reason, 
no  error  analysis  will  be  given  here. 

Most  gap  calibrations  are  based  on  the  well  known  relation  between 
r.ressu re,  P,  shock  velocity,  U,  the  particle  velocity,  u,  and  the 
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and 


U  =  a  +  b  u  (2 ) 

where  a  and  b  are  known  constants  of  the  gap  material.  The  linear 
relation  between  U  and  u  is  one  way  of  describing  the  Hugoniot  of  a 
material.  For  many  materials,  the  linear  relation  holds  over  a  wide 
range,  starting  at  low  pressure.  However  the  Hugoniot  of  PMMA  is 
noticeably  curved  at  low  pressure  as  will  be  discussed  in  a  later 
section.  The  linear  relation  will  be  useful  in  considering  the  errors 
inherent  in  the  methods  discussed  in  the  following. 

If  the  Hugoniot  is  known,  measuring  U  as  a  function  of  x  gives  a 
pressure  calibration.  This  is  the  usual  method  of  calibrating  gap 
tests.  The  dependency  of  shock  velocity  on  distance  is  obtained  by 
differentiating  streak  camera  records  which  give  snock  position  vs 
time.  The  errors  involved  in  this  method  can  be  estimated  by  combining 
Equations  (1)  and  (2)  and  taking  the  total  differential,  AP,  and 
dividing  by  P; 


AP  _  /2U-a  \  AH  +  /  a  \  Aa  +  Ab 
P  “  I  U~a  I  "a  (u-al  a  b  (3) 


Signs  on  the  last  two  terms  have  been  changed  from  negative  to 
positive.  Notice  that  the  coefficients  of  the  first  term  is  greater 
than  1.0.  This  means,  of  course,  that  the  error  in  P  is  greater  than 
that  in  U;  in  this  sense,  the  calculation  of  P  amplifies  the  error  in 
U.  Also,  when  U-^a,  this  amplification  is  more  serious  so  that  the 
accuracy  of  the  calibration  at  low  pressure  is  worse  than  that  at 
high  pressure. 

Shock  velocity,  U,  is  obtained  by  either  a  numerical  or  graphical 
differentiation  of  time-distance  data,  usually  from  streak  camera 
records.  Accuracy  of  these  methods  is  difficult  to  estimate;  it  is 
generally  conceded  that  differentiation  of  numerical  data  is  not  a 
highly  accurate  process.  Our  data  are  probably  somewhat  inaccurate 
also,  we  have  two  sources  of  error  in  our  values  of  U.  We  can 
obtain  some  feel  for  the  problem  by  taking  the  difference  between 
values  of  U  from  two  different  sets  of  data  --  see  Section  3.  This 
difference  is  a  maximum  of  0.12  for  U  a  4.0  mm/usec,  or  about  3%. 
Hence  our  error  in  P  from  this  source  is  1] %  because  the  coefficient 
in  Equation  (3)  is 


(2  x  4.0  -  2 . 5 }  /  (4 . 0  -  2.5)  *  3.7 
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The  range  of  U  is  from  6.0  to  3.0  mm/psec  in  our  calibration.  A  3% 
error  in  U  then  contributes  from  8%  to  21%  to  the  error  in  P  over 
this  range.  This  is  the  chief  fault,  of  this  method.  Accurate  values 
of  U  are  very  difficult  to  obtain  c  d  the  error  in  P,  on  a  relative, 
or  oercentage  basis  is  even  larger  an  that  in  U.  In  the  above  it 
was  assumed  that  U  could  be  measuret  to  the  same  percentaae  error 
over  the  range  6.0  to  3.0.  Had  we  a:  >umed  that  the  absolute  error 
is  0.12  over  the  range,  our  percent  error  would  be  28%  in  n  at  the 
lower  value  of  U. 


Equation  (3)  shows  that  the  error  in  a  is  also  amplified,  while 
that  of  b  is  only  added  in  the  accumulation.  Thr'Se  constants  are 
usually  determined  by  a  least  squares  method.  Values  of  a  and  b, 
and  in  most  cases,  their  respective  quadratic  mean  errors  (ome's) 
are  as  follows  tor  some  representative  data. 


Source  of  data 

High  pressure  data  of  Price 
and  Liddiard^ 


70  mm  camera  data  of  Liddiard 
and  Price4 

Hauver  and  Melani5 


b 


2 . 554  ±0 . 017 


1.618*0.019 


2.554-0.015  1.693*0.016 

2.68  1.61 


For  a  given  set  of  data,  the  constants  are  well  determined  —  their 
QME ‘ s  are  small.  The  variation  in  the  values  of  a  is  about  0.13,  or 
5%  of  the  average  value  of  a.  The  variation  in  b  is  also  about  5%. 
These,  when  added  to  the  above  estimate  for  the  contribution  of  the 
error  in  U,  give  a  range  in  the  variation  of  P  of  from  18%  to  31%  for 
U  varving  from  6  to  3  mm/psec.  These  possible  errors  are  greater  than 
desirable  --  other  methods  should  be  examined  in  the  hope  of  obtaining 
a  more  accurate  calibration. 

An  alternate  method  is  to  measure  the  narticle  velocity  as  a 
function  of  x.  Again,  it  is  assumed  that  the  Hugoniot  is  known.  For 
an  analysis  of  the  error  of  this  method,  it  is  also  assumed  that  the 
U,u  relation  is  linear.  Combining  Equations  (1)  and  (2)  as  before  and 
taking  the  total  differential  gives 


AP 

P 


A  a  ub  Ab 
"a  +  TT  T 


(4) 


where  it  has  been  convenient  to  retain  U  in  place  of  the  quantity 
(a  +  bu) .  Here  none  of  the  coefficients  behave  badly.  That  is, 
none  of  them  tend  to  increase  greatly,  which  is  an  advantage  of  this 
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method.  The  coefficients  are  all  funcfc ions  of  u ,  aa  .it.  it*  of  intorttafe 
to  tabulate  them  over  a  range  of  values  of  u.  Thu  .so  are  given  in 
Table  1,  for  which  U  «  2.576  +  1.602  u.  Note  that  the  first  and  third 
coefficients  increase  with  u.  These  do  not  increase*  oat-.aat rophically , 
however.  Assume  that:  Au  «  0.01,  and  that  Aa/a  »  Ab/b  «  0.05,  or  3%, 
Note  that  over  the  range  0.1  £  u  £  1,8,  the  error  in  u  rnngoa  from 
10%  to  0.6%.  The  errors  in  P  (in  U  are  given  in  Table  l  under  the 
heading  .  These  range  from  16?.  down  to  3.7%  while  U  varies  from 
2.74  to  6.1.  This  is  considerably  bettor  than  the  results  for  the 
preceding  case. 

The  values  assumed  for  the  errors  in  u,  a  and  b  may  not  be 
realistic;  sinal1  values  could  have  been  used  for  him  errors  in  the 
constants,  fo>*  example.  The  values  which  wore  used  are  consistent 
with  those  used  in  the  previous  example,  however.  In  tlu*  case  of  u, 
the  error  may  be  larger  than  that  assumed.  In  Section  4  it  will  be 
shown  that  replicate  shots  give  a  spread  as  groat  as  0,03  and 
0.04  mm/ysec.  These  values  would  increase  the  estimated  error  in  P 
considerably.  We  can  hope  that  tho  curvo  we  draw  through  the  points 
on  a  plot  of  u  vs  x  gives  values  of  u  which  are  more  accurate  than 
0.04  in  absolute  value  over  most  of  the  range.  That  is,  the  smoothing 
which  results  when  a  curve  is  fitted  to  the  points  reduces  tho  error 
in  u  at  individual  values  of  x. 

A  disadvantage  of  measuring  u  vs  x  for  gap  calibration  is  that 
many  experiments  must  be  done  to  establish  tho  relation.  The  relation 
turned  out  to  be  somewhat  more  complicated  than  was  anticipated,  as 
will  be  pointed  out  in  Section  4.  In  contrast,  the  method  based  on 
measuring  U,  requires,  in  principle,  only  a  single  streak  camera 
record  for  the  calibration.  The  work  has  rarely  been  this  easy;  it 
is  found  that  each  record  gives  a  somewhat,  different  curve,  possibly 
due  to  variations  in  the  experiment,  as  well  as  in  reading  and  inter¬ 
preting  the  record.  Usually  the  average  of  several  3hots  is  used  for 
a  calibration.  But  far  fewer  shots  are  required  than  when  particle 
velocity  is  measured. 

Finally,  both  u  and  U  can  be  measured  as  functions  of  x.  Then 
the  error  in  P  is 


\u  AU 
u  U 


The  error  in  P  in  %  is  given  in  the  last  column  of  Table  l,  based  on 
Au/u  =  0.01  and  AU/U  =  0.05  (see  preceding  example).  These  errors, 

E 2  in  Table  1,  are  comparable  to  those  labeled  El .  We  have  the 
alternatives  —  if  the  Hugoniot.  is  well  known,  measuring  u  is  suffi¬ 
cient.  Otherwise,  measure  both  u  and  U  as  functions  of  x.  A  variant 
of  this  latter  method  was  used  when  free  surface  velocity  was  measured 
as  a  function  of  the  gap  length. 4  The  calibration  of  this  report  is 
a  variant  also.  We  measure  both  u  and  U  as  functions  of  x,  extract 
Hugoniot  data  from  the  measurements,  then  use  the  u  vs  x  relation  to 
obtain  P  vs  x. 


nom'r  ‘M-ia 


Table  l 


Krrera 

When  Partiele 

Veioainiea 

Are  1 

Measured 

u 

U* 

(1  ♦  bU/U) 

a/U 

bu/U 

K1 

mm/ w a© o 

n\W/WBOQ 

* 

r 

0.1 

2.74 

1.09 

0,94 

0,09 

16,0 

13,0 

0.2 

2.90 

1.11 

0,09 

0,11 

10.6 

10,0 

0,6 

.1.54 

1,27 

0,73 

0,27 

7,1 

6,7 

1,0 

4.10 

1.10 

0,62 

0,  36 

6,4 

6,0 

1,4 

4.02 

1,47 

0,31 

0,47 

6,0 

3,7 

1.8 

5.46 

1,51 

0.47 

0,33 

3,8 

3.6 

2.0 

5,78 

1.55 

0,43 

0,33 

5.8 

3.3 

2.2 

6.10 

1.50 

0,42 

0,50 

5,7 

5 . 4 

*U  *  2.576  +  1.602  u. 
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NObTR  7.1- 1  ft 

Tho  method  lUscussuri  abuvu  gives  this  Hugonlut  of  the  gap  material, 
Aa  will  In*  shown  in  Heel ion  ft,  the  resulting  Huqoniot  of  PMMA  has 
some  peculiar  fwatursH  which  wo  attribute  in  errors  in  our  moosuro- 
wont n i  Examining  this  relation  has,  hopefully*  helped  ua  tv*  reduce 
some  nf  the  errors,  The  errors  are  pronably  caused,  At  least  in 
part,  by  the  rapid  Attenuation  nf  the  shook  wave,  amt  the  interaction 
of  relief  waves  from  the  periphery  of  the  cylinder,  in  work  with 
the  primary  objective  of  HugoniuV-  measurements,  both  of  these  oompl i - 
cations  are  excluded  from  the  experiments  by  proper  design.  Hence  it 
ia  advantageous  to  use  a  material  in  the  gap  which  has  al ready  been 
characterised  in  careful  studios.  Then  the  preferred  method  for 
calibrating  the  yap  teat  would  be  to  measure  the  particle  velocity  as 
a  function  of  x,  using  the  electromagnetic  pace  as  has  been  done 
in  this  work. 


3.  SHOCK  VELOCITY  IN  PMMA  SHOCKED  BY  PENTOUTE 

This  section  gives  an  account  of  the  measurement  of  the  shook 
velocity  in  the  PMMA  attenuator  of  the  l*SQT  when  pentolite*  is  used 
aa  the  donor  explosive,  in  the  past,  moat  pep  teats  have  boon 
calibrated  by  the  procedure  described  in  the  remainder  of  this 
section,  That  la,  the  shock  velocity  was  determined  as  a  function 
of  the  length  of  the  attenuator)  then  the  llugoniot  equation  of  state 
of  the  attenuating  material  was  uaud  to  convert  shook  velocities 
to  pressures.  Now,  we  are  not  so  dependent  on  the  shock  velocity 
measurements i  independent  measurements  of  narticle  velocity  have 
also  been  made .  Some  of  the  results  of  this  section  will  be  used 
in  conjunction  with  particle  velocity  measurements  to  establish  a 
better  lluqoniot  for  PMMA,  see  Section  5. 

Experimental .  The  experimental  arrangement  for  measuring  shook 
velocity '“In  "tile  PMMA  attenuator  is  given  by  Liddiard  and  Price**  and 
others.  Essentially  the  same  arrangement  was  used  in  this  work. 

The  exception  is  that  no  baffle  was  used  between  the  pentol.it©  and 
the  PMMA,  The  PMMA  was  prepared  bo  that  the  distance  between  flats, 
see  Figure  l,  was  50.8  mm  or  slightly  greater.  This  arrangement 
seems  to  be  adequate,  at  least  in  most  cases,  to  prevent  the  doto- 
nation  product  gases  from  interfering  with  photographing  the  shock 
front . 


*Tha  pentolito  pellets  wore  obtained  from  NAD,  Crane,  Indiana.  They 
replace  graphited  fcetryl  pellets,  «  1.51  g/cc,  wh:  eh  are  no  longer 
available,  The  pentolito  pellets  are  .identified  as  Federal  Stock 
No.  1375-991-8891.  They  are  described  in  BUWEPS  drawing  No,  210  8395, 
"Pellet,  Pentolite,  Donor,  Standard",  and  are  pressed  to  a  density 
of  1.56  g/cc  as  per  "Advisory  Process  For  The  Pressing  Of  Standard 
Donor  Pentolite  Pellets",  NAWEPS  OD  29872,  16  June,  1964.  See 
the  Appendix  for  comments  on  the  uniformity  of  the  dimensions  of 


the  pellets. 
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Span  If icfttionM  for  proparing  tho  .u  t  onunt.ion  cylinders  are 
given  in  Figure  1..  The  cylinders  wore  machined  and  pnliahtul  in  the 
Optics  Shop  of  the  National  bureau  of  Standards  starting  with  63  mm 
(2.5  inch )  diameter  east  rod,  They  wore  free  from  the  optical 
distortion  reported  on  cylinders  uuod  in  earlier  work.  That  is,  the 
flats  were  not  rounded  during  the  polishing  operation  to  form  a 
uarrol  lens,  Also,  the  ends  of  the  flats  were  kept  parallel  to  the 
axis  of  the  cylinders. 

Other  system  improvements  suggested  in  reference  (6)  wore  (a)  use 
of  a  better  lens  in  the  light  saureo,  (b)  glass  window  in  tho  camera 
oort,  (c)  parallel  light  through  tho  sample,  (d)  batter  focusing  of 
the  camera  and  (o)  restricting  the  length  being  observed.  Those 
improvements  will  be  discussed  in  turn. 

Light  for  observing  the  shock  front  in  the  l'MMA  is  obtained  by 
discharging  a  capacitor  through  a  bridgowiro.  The  light  ia  colli¬ 
mated  by  using  a  simple  lens.  In  tho  past,  those  have  beer.  152  mm 
diameter,  356  mm  focal  length  piano  convex  lenses  (tho  same  as  or 
similar  to  Stock  Number  1189  in  tho  Catalog  of  Kdmund  Scientific 
Company),  These  lenses  show  definite  striations  when  observed 
obliquely  in  reflected  light.  Lonsos  used  in  this  work  were  91  nun 
in  diamotor  and  had  a  focal  length  of  215  mm.  These  were  obtained 
from  A.  Jaegers  as  Stock  Number  29A2567,  Although  these  are  not 
high  quality  lenses,  they  do  not  show  the  striations  mentioned  above, 
and  are  therefore  considerably  better  than  those  used  previously. 

Glass  windows  were  used  in  the  camera  port  for  the  shots  reported 
hero.  These  windows  were  cut  from  "Select"  quality  plate  glass 

1.0  inch  thick.  Use  of  glass  windows  gives  demonstrably  better  still 

photographs  in  the  streak  camera.  The  improvement  is  assumed  to 
carry  over  to  the  photographing  of  dynamic  events. 

During  shot  setup,  the  light  source  and  the  collimating  lens 
were  arranged  to  give  parallel  light  through  the  specimen.  The  reason 
for  this  arrangement  was  discussed  in  reference  (6).  The  use  of 
parallel  light  makes  it  impossible  to  photograph  a  10  cm  length  of  the 

specimen  because  the  aperture  in  the  camera  port  is  only  8.8  cm  in 

diameter.  The  field  of  view  was  restricted  to  about  60  mm.  For  two 
shots,  the  section  observed  was  between  0  and  60  mm  from  the  loaded 
end  of  the  specimen.  For  the  second  pair  of  shots,  the  section 
between  about  40  and  110  mm  from  the  loaded  end  was  observed.  As 
noted  in  reference.  (6)  ,  this  arrangement  gives  more  nearly  optimal 
recording  on  the  film.  That  is,  the  slope  of  the  trace  on  the  film 
was  closer  to  unity  than  in  earlier  work.  The  camera  writing  speed 
was  4.0  mm/iisec;  a  speed  of  5.0  mm/ usee  would  have  improved  the 
records  for  the  first  pair  of  shots  (x  60  nun)  . 

The  camera  was  focused  on  the  plane  containing  the  axis  of 
symmetry  of  the  PMMA  cylinder.  The  reason  for  this  is  also  given 
in  reference  (6 )  . 
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Data  Reduction.  The  streak  aamora  records  give  the  position  of  the 
snK>c^”Tn)n'r#  x,  aa  a  function  of  time,  T.  These  records  are  "road'' 
on  a  Universal  Toleroador  ub  described  previously , 6  * 7  There  results 
a  set  of  digitized  data  giving  pairs  of  values  of  x  and  T.  These 
data  wore  smoothed  and  differentiated  using  the  method  described  in 
reference  (6).  This  method  was  called  the  local  smoothing  and  differ¬ 
entiation  (Isd)  method.  The  first  step  of  the  method  consists  of 
getting  a  "smoothed"  value  of  T  for  a  particular  value  of  x  by  least 
square  fitting  a  quadratic  to  the  four  neighboring  points.  The 
quadratic  is  then  solved  for  the  new  value  of  the  time,  Ts.  Each 
point  in  the  set  of  data  is  treated  similarly  to  generate  a  new  set, 

Ts  vs  x. 

The  x , Tg  set  of  data  is  differentiated  by  fitting  a  quadratic 
to  three  points.  This  gives  an  exact  fit  (except  for  round  off 
error)  to  the  points.  The  shock  velocity  is  determined  from  the 
derivative  of  the  quadratic  at  the  second  of  the  three  points  being 
considered.  There  result  tables  of  x,  Tg  and  U“l.  Finally,  the 
values  of  U”1  are  smoothed  in  the  same  manner  that  was  used  for 
smoothing  values  of  time  except  that  a  straight  line,  rather  than  a 
quadratic,  is  used  tc  fit  five  points  instead  of  four. 

In  the  numerical  work,  x  is  used  as  the  independent  variable 
for  reasons  given  in  reference  (6).  Hence  the  derivative  mentioned 
above  is  the  reciprocal  of  the  shock,  velocity,  i.e.,  dT/dx. 

Results  of  the  first  trial  at  reducing  the  data  for  the  four 
shots  are  given  in  Figure  2  and  Tables  2  through  5.  In  the  figure, 
the  shock  velocity,  U,  is  given  as  a  function  of  the  distance,  x. 

In  the  tables,  the  first  and  second  columns  give  T  and  x,  respectively 
in  units  of  "counts"  on  the  Telereader;  in  the  next  two  columns  they 
are  given  in  psec  and  mm  respectively.  The  fifth  column  gives  the 
shock  velocity  in  mra/psec  as  evaluated  by  processes  discussed  above. 
First  divided  differences  in  x  and  T  (here  T  is  used  as  the  indepen¬ 
dent  variable)  are  given  in  the  last  column  of  each  table.  Divided 
differences  are  useful  in  evaluating  the  quality  of  the  data  in  the 
following  discussion. 

For  Shots  865  and  867,  the  camera  recorded  shock  travel  from 
x  =  0.0  to  about  60  mm.  Results  for  Shot  865  appear  to  be  good 
(comparable  with  the  better  shots  in  reference  (6) ,  for  example) . 

This  judgment  is  based  on  observing  the  behavior  of  U  in  Figure  2, 
and  on  the  divided  differences  in  Table  2.  These  latter  generally 
decrease  with  increasing  x  as  they  should.  There  are  "irregularities" 
in  the  differences;  the  sixth  value  is  greater  than  the  preceding 
value  so  that  AU/Ax  >  0.0.  This  would  not  happen  if  each  pair  of 
values  of  x  and  T  lay  on  a  smooth  curve.  (Note:  The  first  entry  in 
Column  6  has  no  meaning).  Continuing  the  examination,  the  13th  is 
greater  than  the  two  preceding  it.  These  irregularities  in  the  data 
did  not  defeat  the  scheme  used  to  obtain  the  values  of  U  —  note  the 
results  in  Column  5,  Table  2  and  in  Figure  2. 
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TABLE  2 


DATA  AND  RESULTS  FOR  SHOT  865 


TIME 

DISTANCE 

TIME 

DISTANCE 

U 

DIFF. 

COUNTS 

COUNTS 

msec 

MM 

MM /mS  EC 

MM/mSEC 

0.0 

16,0 

.000 

.255 

8.783 

0.000 

10.0 

39.0 

.038 

.622 

7,664 

9.769 

19.0 

55.0 

.072 

•  878 

7,041 

7.568 

35,0 

76.0 

.131 

1.213 

6,511 

5.599 

59.0 

106,0 

.221 

1.691 

6.130 

5.334 

90.0 

153.0 

.337 

2,441 

5.880 

6.462 

136.0 

214.0 

.509 

3.414 

5.764 

5.657 

174.0 

267.0 

.651 

4.260 

5.665 

5.948 

219.0 

327.0 

.820 

5.217 

5.521 

5.687 

266.0 

390.0 

1.003 

6.222 

5.379 

5.485 

316.0 

447.0 

1.182 

7.132 

5.292 

5.068 

371.0 

512.0 

1.388 

8.169 

5.212 

5.044 

455.0 

617.0 

1.702 

9.844 

5.127 

5.334 

546.0 

725.0 

2.042 

11.567 

5.031 

5.066 

643.0 

844.0 

2.405 

13.466 

4.920 

5.235 

771.0 

988.0 

2.883 

15.763 

4.758 

4.803 

905.0 

1128.0 

3.384 

17.997 

4.617 

4.462 

999.0 

1233.0 

3.735 

19.672 

4.529 

4.769 

1122.0 

1362.0 

4.195 

21.730 

4.445 

4.479 

1309.0 

1549.0 

4.893 

24.714 

4.354 

4.272 

1459.0 

1703.0 

5.453 

27.171 

4.284 

4.385 

1626.0 

1870.0 

6.077 

29.835 

4.208 

4.272 

1812.0 

2052.0 

6.772 

32.739 

4.119 

4.180 

1961.0 

2195.0 

7.328 

35.021 

4.044 

4.100 

2131.0 

2356.0 

7.963 

37.589 

3.949 

4.046 

2276.0 

2486.0 

8.504 

39.664 

3.860 

3.832 

2411.0 

2609.0 

9.008 

41.626 

3.772 

3.893 

2586.0 

2764,0 

9.662 

44.099 

3.651 

3.785 

2715.0 

2872.0 

10.143 

45.8  2 

3.571 

3.579 

2854.0 

2984.0 

10.662 

4  7  ,  b  0  9 

3.499 

3.445 

2974.0 

3084.0 

11.110 

49.205 

3.441 

3.562 

3103.0 

3184.0 

11.591 

50.800 

3.401 

3.315 

3123.0 

3197.0 

11.666 

51.007 

3.394 

2.781 

3222.0 

3284.0 

12.035 

52.395 

3,371 

3.756 

33 63 .0 

3386.0 

12.561 

54,023 

3.350 

3.094 

3479.0 

3480.0 

12.994 

55.523 

3.331 

3.464 

3600.0 

3574.0 

13.446 

57,022 

3.312 

3.322 

3715.0 

3662.0 

13.875 

58.426 

3.289 

3.272 

3790.0 

3720.0 

14.155 

59.352 

3.274 

3.307 

3867.0 

3779.0 

14.442 

60.293 

3.258 

3.276 
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ghot  867  gave  results  which  are  inferior  to  those  of  865,  see 
Table  3  and  Figure  2.  Values  of  U  are  too  small  for  x  <  5  mm,  and 
for  x  >  47  nun.  This  "badness"  is  in  the  digital  data  and  presumably 
in  the  streak  camera  trace  itself  —  it  is  not  an  artifact  of  the 
numerical  operations.  This  conclusion  is  substantiated  by  the 
behavior  of  the  differences  in  Table  3.  Values  of  the  differences 
do  not  decrease  with  x  smoothly.  The  same  situation  is  observed 
between  49.6  and  51.3  mm,  and  near  the  end  of  the  table. 

For  Shot  866  (x  varies  from  47  to  109  mm) ,  U  increases  with  x 
for  47  <.  x  <.  51  and  for  75  <_  x  <.  82,  see  Table  4.  The  first  3  values 
of  the  differences  bear  out  the  first  observation.  The  increase  of 
U  with  x  near  80  mm  may  be  due  to  the  data  which  resulted  in  the 
value  3.211  for  the  difference  at  80.3  mm.  For  this;  shot  also,  there 
are  irregularities  in  the  data  which  are  not  due  to  the  numerical 
methods.  Unfortunately,  the  numerical  methods  cannot  remove  the 
irregularities . 

Results  for  Shot  869  (Table  5)  are  much  the  same  as  those  of 
Shot  866  discussed  immediately  above.  There  must  be  a  reproducible 
cause  for  the  peculiar  behavior  observed  at  the  beginning  of  each  of 

these  records.  It  was  thought  that  this  behavior  was  due  to  inaccu¬ 

rate  flats  near  the  ends  of  the  PMMA  cylinders.  More  care  in 
preparing  these  flats  has  not  eliminated  the  trouble.  In  the 
following,  these  irregular  results  (at  the  beginning  of  the  records, 
and  the  oscillatory  behavior  later  in  the  records  of  866  and  868, 
see  Figure  2)  are  treated  as  errors,  and  the  deviations  from  the 
expected  results  are  eliminated,  or  at  least,  smoothed. 

Defects  in  the  results  could  possibly  be  reduced  by  using  more 
forceful  smoothing  of  the  x,T  and  the  x,U-1  data  using  graphical  or 
numerical  methods.  It  was  done  by  removing  those  data  which  appear 

to  be  the  chief  causes  of  the  trouble.  This  is  a  subjective  process; 

it  was  guided  by  examining  the  behavior  of  U  and  the  divided  differ¬ 
ences  in  Figure  2  and  in  Tables  2  through  5.  Data  eliminated  are  the 
following ; 


Shot  No .  Data  Removed  For  Values  of  x  of 

865  none 

867  0.0  through  4.54  mm 

49.6  through  52.1  mm 

55.7  through  57.7  mm 

866  47.1  through  49.9  mm 

868  44.2  through  47.9  mm 


The  revised  sets  of  data  were  treated  in  the  same  way  as  the 
originals,  see  above.  Results,  which  are  shown  in  Figure  3  are  less 
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TABLE  3 

DATA  AND  RESULTS  FOR  SHOT  867 


TIME 

DISTANCE 

TIME 

DISTANCE 

U 

DIFF. 

COUNTS 

COUNTS 

,uSEC 

MM 

MM//ASEC 

MM/>SEC 

0.0 

0.0 

0.000 

0.000 

5.618 

0.000 

23.0 

43.0 

.086 

.680 

5.610 

7.912 

66.0 

88.0 

.246 

1.392 

5,601 

4.439 

107.0 

145.0 

.  399 

2.294 

5.596 

5.892 

155.0 

209.0 

.579 

3.307 

5.573 

5.651 

214.0 

287.0 

.799 

4.541 

5.527 

5*604 

276.0 

368,0 

1.030 

5.822 

5.448 

5.538 

330.0 

440.0 

1.232 

6,961 

5.358 

5.651 

405,0 

531.0 

1.512 

8.401 

5.219 

5.144 

488.0 

632.0 

1.821 

9.999 

5.079 

5.159 

551.0 

707.0 

2.056 

11.185 

4.984 

5.048 

658.0 

828.0 

2.456 

13.099 

4.855 

4.796 

733.0 

916.0 

2.735 

14.492 

4.773 

4.975 

882.0 

1080.0 

3.291 

17.086 

4.641 

4.668 

1045.0 

1256.0 

3.899 

19.871 

4.519 

4.580 

1205.0 

1419.0 

4.496 

22.449 

4.420 

4.322 

1335.0 

1556.0 

4.981 

24.617 

4.345 

4.470 

1446.0 

1670.0 

5.395 

26.420 

4.290 

4.357 

1564.0 

1788.0 

5.835 

28.287 

4.234 

4.242 

1682.0 

1904.0 

6.275 

30.122 

4.168 

4.171 

1799.0 

2018.0 

6.711 

31.926 

4.091 

4.134 

1905.0 

2122.0 

7.106 

33.571 

4.010 

4.163 

2014.0 

2227.0 

7.513 

35.233 

3.931 

4.087 

2153.0 

2349.0 

8.031 

37.163 

3.842 

3.725 

2309.0 

2487.0 

8.612 

39.346 

3.755 

3.754 

2457.0 

2617.0 

9,164 

41.403 

3.685 

3.728 

2571.0 

2717,0 

9.589 

42.985 

3.634 

3.723 

2666.0 

2796.0 

9.943 

44.234 

3.587 

3.530 

2774.0 

2889.0 

10.346 

45.706 

3.525 

3.655 

2913.0 

3003.0 

10.864 

47.500 

3.466 

3.481 

3014,0 

3084.0 

11.240 

48.791 

3.461 

3.404 

3083.0 

3137.0 

11.497 

49.629 

3.444 

3.261 

3126.0 

3173.0 

11.658 

50.199 

3.450 

3.553 

3162,0 

3207.0 

11.792 

50.737 

3.398 

4.007 

3174.0 

3211.0 

11.837 

50.800 

3.421 

1.417 

3207,0 

3240.0 

1  1.960 

51.259 

3.367 

3.730 

3277.0 

3292.0 

12,220 

52.081 

3.335 

3.154 

3343.0 

3345.0 

12,466 

52.920 

3.302 

3.409 

3415.0 

3403.0 

12,735 

53.838 

3.307 

3.419 

3499,0 

3465.0 

13,048 

54.818 

3.279 

3.134 

3566.0 

3519.0 

13.297 

55.673 

3.243 

3.421 

3625,0 

3564.0 

13.517 

56.385 

3.163 

3.238 

3688,0 

3612.0 

1  3.752 

57.144 

3.083 

3.235 

3742,0 

3649,0 

13.953 

57.729 

3.023 

2.909 
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TABLE  4 


data  and  results  FOR  shot 


TIME 

COUNTS 


0.0 

65.0 

153.0 

224.0 

311.0 

307.0 

474.0 

588.0 

724.0 

888.0 

1028.0 

1160.0 

1362.0 

1536.0 

1725.0 

1862.0 

2001.0 

2133.0 

2292.0 

2453.0 

2634.0 

2774.0 

2933.0 

3100.0 

3262.0 

3410.0 

3548.0 

3716.0 

38e6.0 

4034.0 

4207.0 

4337.0 

4477.0 

4608.0 

4761.0 

4946.0 

5091,0 

5212.0 

5317.0 


DISTANCE 

COUNTS 


0.0 
50.0 
119.0 
177.0 
248.0 
307.0 
375.0 
464.0 
572.0 
695 . 0 
805.0 
906.0 
1062.0 
1192.0 
1333.0 
1438.0 
1535.0 
1632.0 
1747.0 
1863.0 
1992.0 
2090,0 
2214.0 
2337.0 
2453.0 
2560.0 
2661.0 
2781.0 
2902.0 
3011.0 
3130.0 
3223,0 
3325.0 
3418.0 
3524.0 
3654.0 
3757,0 
3842.0 
3915.0 


TIME 

msec 

0.000 
.24  2 
.571 
.836 
1.160 
1 . 444 
1 . 768 
2.194 
2.701 
3.313 
3,835 
4.327 
5.081 
5.730 
6.435 
6.946 
7.464 
7.956 
8.549 
9.150 
9.825 
10.347 
10.94 0 
11.563 
12.167 
12. 719 
13.234 
13.860 

14.494 

15,046 

15.692 

16.176 

16.699 

17.187 

17.758 

18.448 

18.988 

19.440 

19.831 


DISTANCE 

MM 

47.130 

47,921 

49.013 

49.930 
51.054 

51.987 
53.063 
54.471 
56.179 
58.125 
59.866 
61.463 

63.931 

65.988 
68.219 
69.880 
71.415 
72.949 
74, 769 
76.604 
78.645 
80.322 
82.157 
84.103 
85.938 
87.631 
89.229 
91.127 
93.041 
94,766 
96,649 
98.120 
99. 734 

101.205 
102.882 
1 04.939 
106.568 
107.913 
109,068 


866 

U 

MM/mSEC 

3.291 

3.314 

3.346 

3.362 

3.364 

3.353 

3.334 

3.312 

3.293 

3.277 

3.262 

3.246 

3.216 

3.184 

3.145 

3.116 

3.094 

3.080 

3.077 

3.082 

3.092 

3.097 

3.097 

3.088 

3.077 

3.066 

3.056 

3.045 

3.036 

3.029 

3.022 

3.017 

3.010 

3.004 

2.996 

2.98  3 

2.980 

2.978 

2.977 


DIFF. 

MM/juSE'C 


0.000 

3.262 

3.323 

3 , 464 

3.461 

3.292 

3,315 

3.*11 

3.368 

3,iei 

3,332 

3.245 

3.275 

3  o  169 

3.164 

3.251 

2.961 

3.117 

3.068 

3.056 

3.023 

3.211 

3.095 

3.124 

3.038 

3.067 

3.105 

3.030 

3.019 

3.124 

2*918 

3.035 

3.091 

3.012 

2.939 

2.981 

3.013 

2.980 

2.950 
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TABLE  5 

DATA  ANO  RESULTS  FOR  SHOT  868 


time  distance  time 

COUNTS  COUNTS  fJ-SEC 

0.0  0.0  0.000 

76.0  58.0  .283 

139.0  113.0  ,518 

211.0  172.0  .786 

273,0  226.0  1,017 

343.0  279,0  1.278 

403.0  326.0  1.501 

469.0  382.0  1,747 

529.0  427.0  1.971 

594.0  476.0  2.213 

689.0  552.0  2.567 

792.0  632 . 0  2.95© 

896.0  712.0  3.338 

996.0  789.0  3.710 

1111.0  875.'  4.139 

1257.0  993.0  4.682 

1389,0  1092.0  5,174 

1541.0  1204.0  5.740 

1707.0  1331.0  6.359 

1868.0  1453,0  6.958 

2n 1 3,0  1559,0  7,498 

2155.0  1664.0  8,027 

2282.0  1756.0  8.500 

2402.0  1843.0  8.947 

2519.0  1928.0  9.383 

2647.0  2021,0  9.860 

2816,0  2146,0  10.489 

2996.0  2282.0  11.160 

3138.0  2386.0  11,689 

3326.0  2522.0  12,389 

3457.0  2618.0  12.877 

3580,0  2707,0  13.335 

3688,0  2783.0  13,737 

3850.0  2901.0  14.341 

3948.0  2973.0  14,706 

4088,0  3072.0  15.227 

4205,0  3156,0  15.663 

4377.0  3273,0  16.303 

4528.0  3380.0  16,866 

4668,0  3479,0  17.387 

4775,0  3556,  C  17, ,86 

4921,0  3663.0  18.329 

5041.0  7746,0  18,776 

5174,0  3839.0  1 9,272 

5267,0  3907,0  19,618 

5323.0  3944,0  19,827 

5380,0  3985,0  20.039 


DISTANCE  U  DIFF, 

MM  MM/^tSEC  MM/jixSEC 

44.248  3.459  0.000 

45.172  3.475  3.266 

46.049  3.489  3.735 

46.989  3.489  3.506 

47,850  3.470  3.726 

48.695  3.436  3.240 

49.444  3.404  3.352 

50.337  3.371  3.630 

51.054  3.357  3.210 

51*835  3.341  3.226 

53.046  3.325  3.423 

54.322  3.312  3.324 

55®597  3.308  3.292 

56.824  3,298  3.295 

58.195  3.288  3.200 

60.076  3.271  3.458 

81.654  3.257  3.210 

63.439  3.235  3.153 

65.464  3.208  3.274 

67.408  3.181  3.243 

69.098  3.158  3.128 

70.771  3.137  3.164 

72*23 9  3.124  3.100 

73.625  3.122  3.103 

74,980  3.129  3.109 

76.462  3.138  3.109 

78.454  3,147  3.165 

80.622  3.145  3.233 

82.280  3.136  3.134 

84.448  3.115  3.096 

85.978  3.102  3.136 

87,397  3.094  3.097 

88.608  3.087  3.012 

90.489  3.074  3.117 

91.637  3.059  3.144 

9  3 . 2 1  ¥>  3,040  3.026 

*’4.554  3.028  3.073 

96,419  3.024  2 . ° 1 1 

98.124  3.032  3.033 

99,702  3.041  3.026 

100.92°  3,04fc  3.080 

102,635  3.042  3,136 

103.958  3,033  2.960 

105.440  3.018  2.993 

106.524  3,015  3.129 

107.114  3.014  2.828 

107,768  3.012  3.078 


DISTANCE  ,  mm 

FIG .  3  SHOCK  VELOCITY  VS  DISTANCE,  FINAL  RESULTS 
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erratic  than  those  of  Figure  2,  as  was  expected.  It  is  now  necessary 
to  combine  all  the  data  to  give  a  single  U  vs  x  curve.  This  could 
be  done  graphically.  Instead,  a  new  set  of  results  giving  U  at  each 
mm  in  x  for  each  shot  are  obtained  by  interpolating  in  the  sets 
described  above.  Representative  samples  of  these  new  sets  are  given 
in  Table  6.  The  last  column  gives  the  differences  in  U  for  replicate 
shots;  these  tell  something  about  the  precision  of  the  work.  Note 
that  these  differences  are  greatest  at  35  and  40  mm.  This  is  the 
middle  region  of  the  "hump"  previously  observed  in  the  U  vs  x  curve. ^ 
It  may  be  caused  by  interaction  of  rarefaction  waves  with  the  shock 
wave  as  has  been  mentioned  above.  This  interaction  apparently 
introduces  some  imprecision  in  our  measurement  of  shock  velocity. 

Having  values  of  U  at  equally  spaced  values  of  x  for  each  pair 
of  shots  permits  us  to  combine  the  results  by  averaging.  Results 
for  Shots  865  and  867  are  averaged  for  6  <.  x  49;  Shots  866  and  868 
are  averaged  for  50  <.  x  <1  105.  These  new  results  are  displayed  in 
Figure  4. 

The  oscillations  in  the  U,x  curve  of  Figure  4  may  or  may  not  be 
real  —  they  are  in  the  digital  data,  an  observation  borne  out  by 
the  behavior  of  the  divided  differences.  Here  it  is  assumed  that 
AU/Ax  should  be  negative  for  all  values  of  x.  This  condition  can  be 
forced  on  the  results  by  fitting  the  averaged  data  (see  above)  with 
a  cubic  for  40  £  x  £  105  mm.  There  results 

U  =  7.416  -  0.151  x  +  0.00179  x2  -  0.00000713  x3,  (6) 


which  fits  the  data  well,  see  the  solid  curve  in  Figure  4.  The 
cubic  has  an  inflection  point  (second  derivative  is  zero)  at  83.5  mm; 
it  appears  to  cause  no  trouble.  As  desired,  the  first  derivative 
remains  negative  over  the  interval  of  x  in  the  data.  For  the  final 
set  of  U,x  data,  the  cubic  is  used  for  45  <.  x  <.  105  mm.  For 
8  <  x  i  45  mm  the  averaged  data  are  used  --  these  data  are  not  fitted 
to  a  function.  The  use  of  these  intervals  insures  that  the  two  sets 
of  data  join  smoothly.  These  data  are  all  shown  in  Figure  4;  a 
sampling  is  shown  in  Column  2  of  Table  7. 

If  we  relied  on  measurements  of  shock  velocities  for  the 
calibration  of  the  gap  test  for  pentolite,  we  would  have  to  obtain 
results  in  the  interval  0.0  <.  x  <  8  mm.  Because  the  shock  velocity 
decreases  with  x  very  rapidly  in  this  interval,  it  is  much  more 
difficult  to  measure  than  elsewhere.  In  the  past,  the  Chapman-Jouguet 
variables  have  been  evaluated  for  the  explosive.  Then  a  shock  imped¬ 
ance  method  was  used  to  calculate  the  pressure,  particle  velocity 
and  shock  velocity  at  the  interface.  This  procedure  ignores  any 
effect  of  the  finite  reaction  zone  in  the  explosive.  Here,  we  do 
not  really  need  the  shock  velocities  close  to  the  explosive  because 
we  have  particle  velocity  measurements  instead  (see  next  Section) . 
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Table  6 

SPREAD  IN  SHOCK  VELOCITY  BETWEEN  SHOTS 
(Interpolated  Data) 


_ l?i  mm/ usee _ 

x  Shot  Shot'  Shot  shot;"  Difference 

mm  865  867  866  868 


5 

5.551 

10 

5.124 

5.065 

15 

4.807 

4.732 

20 

4.498 

4.504 

25 

4.346 

4.333 

30 

4.208 

4.173 

35 

4.045 

3.941 

40 

3.851 

3.730 

45 

3.604 

3.562 

50 

55 

60 

65 

70 

75 

80 

85 

90 

95 

100 

3.428 

3.388 

0.059 

0.075 

0.006 

0.013 

0.035 

0.104 

0.120 

0.042 

0.040 


3.300 

3.313 

0.013 

3.264 

3.272 

0.008 

3.201 

3.214 

0.013 

3.108 

3.143 

0.035 

3.073 

3.130 

0.057 

3.09  7 

3.148 

0.051 

3.083 

3.107 

0.024 

3.050 

3.080 

0.030 

3.028 

3.016 

0.012 

3.010 

3.048 

0.038 

18 


NOLTR  n-u 


Table  7 

COMPARISON  01'  SHOOK  VELOCITIES  PROM  TURKS 
SETS  OP  SHOTS  USING  PENTOLITE 
DONOR  AND  ONE  SET  USING  TKTRYL  DONORS 


Distance 

mm 


mr 

work 
(Lot  3) 


Shock  Velocity 
mm  At  set? 


And  716 
{tot  2) 


Through  380 
(tot  2) 


T$ 

Calibra- 
tion  m 


Calibra¬ 
tion  U) 


0,0 

6.24 

6.35 

5.9  5 

(6.00) 

(6.24) 

3,0 

5,33 

3.41 

5.34 

(5.39) 

(5.58) 

10.0 

5.10 

5.05 

4.93 

4.94 

5.09 

13.0 

4,77 

4.78 

4  .60 

4.63 

4.76 

20.0 

4,30 

4.37 

4.43 

4.39 

4.46 

25.0 

4,34 

4,33 

4  . 2  3 

4.19 

4.22 

30,0 

4.19 

4.18 

4.11 

4.01 

4.04 

35,0 

3.99 

3.96 

3.90 

3.04 

3.84 

40,0 

3,77 

3.77 

3.71 

3.66 

3.66 

45.0 

3.50 

3.58 

3.58 

3.50 

3.50 

50,0 

3.43 

J  .  4U 

J  .  •*  2 

3.40 

3 . 40 

55.0 

3.32 

3. 28 

3.30 

3.34 

3.34 

60.0 

3.23 

3.23 

3.27 

3.28 

3.28 

65.0 

3.17 

3.22 

3.26 

70,0 

3.14 

3.15 

3.2  5 

3.20 

3.20 

75,0 

3.11 

3.13 

3.21 

80.0 

3.10 

3.11 

3.16 

3.15 

3.15 

83.0 

3.09 

3,13 

3.16 

90,0 

3.08 

3.12 

3.15 

3.12 

3.12 

95,0 

3.06 

3.02 

3.14 

100.0 

3.02 

2,98 

3.06 

3.10 

3.10 

NOTE  t 

Second  column 

is  for  our 

new  stock 

of  pentolite. 

Third  and  fourth  columns  arc  for  pentolite  described  in 
reference  (1)  . 


Shot  715  -  record  covered  0  t  x  1  55  mm 
Shot  716  -  covered  ~<5  i  x  i  105  mm 

Shota  578  through  580  -  records  covered  0  s.  x  ±  110  mm 
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Comparison  With  Proviouu  Rnniiltn.  At  thlw  point,  it  ia  of  Intoroafc 
to  uomparo  tins  roauiln  from  Uv  piouuttinq  auction  with  Home  of  those 
obtained  previously.  The  new  uault.s  are  gi  von  in  Table  7,  Column  a  1 
and  a»  well  an  in  Figure  <4  Those  results  in  Columns  3  through  6 
will  bo  discussed  in  turn  in  .a,  fnllowinq.  Thu  explosive  used  in 
the  shots  represented  by  the  tt.*  in  Columns  3  and  4  is  pentolite 
from  Lot  2,  now  references  (1.)  a*  (o). 

Column  3  in  Tahiti  7  gives  the  results  from  tho  "close-up"  shots 
described  in  reference  (6),  For  wit  713,  tho  camera  viewed  abou- 
50  mm  of  tho  specimen  adjacent  to  tie  explosive.  For  Shot  716,  the 
camera  viowod  tho  region  from  about  30  to  100  mm  from  tho  explosive. 
Values  of  U  in  Column  3  wore  obtain  id  by  smoothing  and  differen¬ 
tiating  tho  data  from  tho  records  an  described  above.  The  discrep¬ 
ancies  between  values  of  1)  in  Columns  2  and  3  and  Columns  2  and  4 
for  \  i  10  mm  are  of  no  consequence.  The  reason  for  this  is  that 
no  attempt  was  made  to  force  the  results  given  in  Columns  3  and  4 
to  tho  selected  value  of  U  at  x  »  0.0.  Agreement  over  most  of  the 
range  of  x  is  good.  Note,  however,  the  value  of  U  at  85  mm  in 
Column  3  is  greater  than  the  value  at  00  nun.  Such  "oscillations" 
wore  removed  from  the  results  in  Column  2  by  fitting  a  cubic  to 
the  data  for  x  greater  than  about  44  mm.  The  two  columns  are  in 
agreement  to  *1%  for  almost  all  values  of  x  >  10  mm.  This  means  that 
the  two  lots  of  pentolite  are  indistinguishable. 

The  values  of  U  in  the  fourth  column  of  Table  7  are  from 
Shots  578,  579  and  580,  They  were  called  "regular  shots"  in 
reference  (G);  for  these  the  camera  recorded  a  field  of  view  somewhat 
qreater  than  100  mm.  These  values  agree  with  those  in  Column  2 
less  well  than  do  those  of  Column  3.  Note  that  these  regular  shots 
give  values  of  U  which  are  greater  than  those  in  the  preceding  two 
columns  for  x  >  55  mm.  Pentolite  from  Lot  2  was  used  for  the  shots 
which  yielded  the  data  in  Columns  3  and  4.  If  we  relied  on  the 
results  from  the  regular  shots  (Column  4)  we  would  conclude  that 
pentolite  from  Lot  2  was  more  energetic  than  that  from  Lot  3, 

(Column  2).  It  is  probable  that  the  data  from  the  close  up  shots 
(Column  3)  are  the  more  accurate,  and,  as  noted  above,  the  two  lots 
of  pentolite  are  indistinguishable  when  the  data  are  obtained  and 
reduced  in  the  same  manner. 

The  numbers  in  Column  5  of  Table  7  represent  the  last  calibration 
of  the  LSGT  with  tatryl.  For  x  >.  55  mm,  the  shock  velocities  of 
tetryl  are  greater  than  those  given  for  the  new  lot  of  pentolite  in 
Column  2.  Because  pentolite  is  the  more  powerful  explosive,  the 
opposite  is  expected.  Again,  it  is  probable  that  the  new  data  are 
the  more  accurate?  if  so,  the  earlier  tetryl  calibration  curve  gives 
values  of  U  which  are  too  large  for  x  >,  55  mm.  This  difference  was 
one  of  the  reasons  for  using  the  electromagnetic  gage  to  measure 
particle  velocities  in  PMMA  attenuators  when  shocked  by  pentolite 
and  by  tetryl. 
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Comparison  of  Column  2  with  Column  5  shows  that  our  present 
values  of  U  for  pentolite  are  greater  than  those  for  tetryl  for 
x  <  50  nun .  This  indicates,  of  course,  that  pentolite  (1.56  g/cc)  is 
more  powerful  than  tetryl  (1.51  g/cc),  and  that  the  P  vs  x  curve  of 
pentolite  will  lie  above  that  of  tetryl.  The  above  remarks  are 
based  on  the  assumption  that  each  set  of  results  is  accurate. 

Because  different  techniques  were  used  in  obtaining  the  two  sets  of 
data,  one  set  may  be  more  accurate  than  the  other.  Again,  measure¬ 
ments  of  particle  velocity  in  the  PMMA  attenuator  help  resolve  the 
issue. 

Because  particle  velocities  have  been  measured  in  the  PMMA 
attenuator  for  both  explosives,  there  is  no  point  in  converting  the 
values  of  U  obtained  in  this  section  to  pressures.  These  data  will 
be  used  in  a  later  section  to  help  establish  a  Hugoniot  for  PMMA, 
which,  along  with  the  particle  velocities,  will  determine  the  pressure 
as  a  function  of  distance. 


4 .  PARTICLE  VELOCITY  MEASUREMENTS 

Particle  velocity,  u,  in  the  PMMA  attenuator  is  measured 
directly  using  the  electromagnetic  velocity  ( EMV )  gage.  The  gage 
consists  of  an  aluminum  foil  (0.0005  inch  thick)  as  shown  in 
Figure  5a.  The  gage  is  formed  by  wrapping  the  foil  around  piece  A 
where  the  dimension  l  is  carefully  measured.  The  part  of  the  foil 
around  the  bottom  of  piece  A  is  called  the  base.  Pieces  B  and  C 
are  cemented  (using  chloroform)  to  piece  A  and  the  foil  and  to  a 
PMMA  cylinder  whose  length  has  been  measured.  This  whole  assembly 
is  then  placed  under  2000  psi  for  20  minutes  to  remove  any  air 
bubbles.  The  booster  pellets  and  the  detonator  holder  are  then  added 
to  the  assembly  which  is  placed  in  a  magnetic  field.  The  assembly 
is  oriented  as  shown  in  Figure  5a,  so  that  the  shock  from  the  donor 
charge  causes  the  base  of  the  gage  to  move  across  lines  of  magnetic 
flux,  generating  an  emf 


V  =  H-u-Jl-10"4  (7) 


where  the  units  are  volts,  gauss,  mm/usec  and  mm,  respectively. 

Note  that  the  legs  of  the  gage  do  not  generate  an  emf.  The  legs 
of  the  foil  are  attached  to  a  coaxial  cable,  see  Figure  5b.  The 
50  ohm  resistor  is  used  to  impedance  match  the  gage  to  the  cable; 
the  oscilloscope  end  of  the  cable  is  also  impedance  matched  through 
a  voltage  divider  so  that  two  oscilloscopes  can  be  used. 

Typical  oscilloscope  records  are  shown  in  Figures  6a  through  6c. 
For  the  record  in  Figure  6a,  the  PMMA  cylinder  was  5  mm  long.  The 
rise  time  for  this  record  was  very  short,  approximately  20  nano¬ 
seconds  (ns).  After  350  ns,  the  record  shows  some  noise.  This  is 
probably  induced  by  breaking  of  the  PMMA  or  by  the  explosive  gases 
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expanding  around  the  edges  of  the  PMMA  cylinder.  A  cardboard  gas 
|  shield  fitted  to  the  PMMA  cylinder  retards  these  gases  sufficiently 

so  that  the  noise  does  not  destroy  the  usefulness  of  the  early  part 
jr  of  the  record. 


m- 


i:: 


The  record  shown  in  Figure  6b  is  from  a  shot  in  which  the  PMMA 
cylinder  was  35.2  mm  long.  The  rise  time  is  about  40  ns.  The 
curious  feature  of  the  record  is  the  overshoot  which  lasts  about 
40  ns.  It  is  doubtful  if  this  is  a  function  of  the  circuitry.  It 
may  be  due  to  the  behavior  of  PMMA  under  rapid  loading  —  the  material 
may  show  relaxation  effects.  It  is  possible,  of  course,  that  wave 
interactions  actually  produce  such  a  velocity  profile.  The  over¬ 
shoot  has  been  ignored  in  the  determination  of  the  peak,  particle 
velocity. 

Figure  6c  is  a  2  microsecond  record  of  the  particle  velocity 
for  a  100.1  mm  long  cylinder.  Rise  time  is  again  short,  and  there 
is  a  slight  overshoot.  The  remainder  of  the  record  indicates  that 
the  particle  velocity  is  relatively  constant. 


£ 


Records  are  read  and  digitized  by  using  a  Universal  Telereader, 
after  which  plots  are  made;  see  Figure  7  for  plots  of  the  records 
shown  in  Figures  6a  -  6c.  The  peak  particle  velocities  can  be 
inferred  by  a  linear  extrapolation  of  the  curves  to  zero  time.  As 
noted  above,  the  overshoot  on  a  record  is  ignored.  Usually,  the 
extrapolation  is  based  on  a  visually  selected  section  of  the  record 
about  0.5  microseconds  long.  Most  of  the  results  were  obtained 
by  then  fitting  this  section  by  least  squares  to  a  straight  line 
over  the  selected  interval  of  the  record. 

Peak  particle  velocities  for  tetryl  donors’^  are  given  in  Table  8 
along  with  the  interval  of  time  of  the  record  used  in  obtaining  the 
oeak  velocity.  Also  included  in  the  table  are  comments  concerning 
some  of  the  records.  Asterisks  following  the  shot  number  means  that 
the  record  is  from  a  Hewlett-Packard  oscilloscope.  The  time  base 
for  these  records  is  usually  5  microseconds  so  that  they  show  less 
'v-trly  detail  than  those  from  the  Tektronix  oscilloscope,  for  which 
the  time  base  is  usually  1  microsecond. 

All  records  obtained  are  listed  in  Table  8,  but  some  of  these 
records  were  unacceptable  because  of  noise  or  form  (116,  121,  244, 
241*)  and  were  therefore  discarded.  Other  records  which  were  not 
used  fall  into  two  groups:  (1)  those  obtained  at  0  <_  x  <  0.86  mm 
and  (2)  those  obtained  with  the  Hewlett-Packard  (HP175)  oscilloscope. 
We  found  in  earlier  work  that  particle  velocities  above  the  theo¬ 
retical  (2.20  mm/ nsec  at  x  =  0  for  tetryl  loading  of  PMMA)  were 
measured  at  small  or  no  attenuation,  and  attributed  this  to  the 
effect  of  the  reaction  zone  in  tetryl. ®  The  reaction  zone  was  not 
considered  in  deriving  the  theoretical  value;  instead  the  C-J 


i'Graphited  tetryl  pressed  to  a  density  of  1,51  ±  0.1  g/cc;  it 
contains  0.5%C. 
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T ABLE  8 

PEAK  PARTICLE  VELOCITIES  FOR  TETRYL  DONOR 


SHOT 

X 

u 

TIME  INTERVAL  COMMENTS 

NO 

MM 

MM/>aSEC 

FOR  LS  FIT 

289 

0.0 

2.11 

0.05-0.16 

POINT  INITIATED 

185 

0.0 

2.39 

1  MIL  FOIL,  PLANE  WAVE  BOOSTER, 

25  MM  OF  TETRYL 

106 

0.25 

2.31 

SOME  NOISE  ON  RECORD,  NOT  USED 

118 

0.25 

2.26 

NOT  USED,  RECORD  ONLY  0.2  USEC 

LONG 

120 

0.25 

2.21 

NOT  USED 

108 

0.86 

2.25 

183 

0.86 

2.22 

1  19 

2.0 

2.05 

117 

3.0 

1.88 

0.02-0.22 

116 

4.0 

1.61 

0.03-0.22 

NOISY  RECORD ,  NOT  USED 

241 

5.0 

1.74 

0.05-0.34 

NOISY  RECORD,  NO  GAS  SHIELD 

244 

5.0 

1.78 

0.11-0.59 

NOISY  RECORD,  NO  GAS  SHIELD 

241# 

5.0 

1.76 

0.35-1.10 

NOISY  RECORD,  NOT  USED 

258 

5.0 

1.75 

0.04-0.34 

GAS  SHIELD  REDUCED  NOISE 

2  42 

10.0 

1  .49 

0.15-0.64 

NOISY  RECORD 

240 

15.0 

1.30 

0.16-0.33 

0.2  USEC  RISE  TIME 

240* 

15.0 

1.38 

0.40-1.50 

NOT  USED,  COULD  ALSO  BE  READ  AS 

1.30  ON  PLATEAU 

121 

20.0 

l  .14 

0.02-0.20 

SHORT  RECORD,  FLAT,  NOT  USED 

243 

20.0 

1.19 

0.13-0.64 

239* 

25.0 

1.11 

0.05-0.46 

OVERSHOOT 

231 

30.0 

1.01 

0.06-0.56 

OVERSHOOT 

231* 

30.0 

1  .03 

0.30-1 .00 

OVERSHOOT,  NOT  USED 

222 

35.0 

0.892 

0.22-0.72 

222* 

35.0 

0.912 

0.20-0.80 

NOT  USED 

230 

35,2 

0.904 

0.08-0.60 

OVERSHOOT 

230* 

35.2 

0.91 

0.08-0.90 

OVERSHOOT,  NOT  USED 

232 

39.9 

0.735 

0.05-0.55 

232* 

39.9 

0.741 

0.05-0.90 

NOT  USED 

218 

45.0 

0.625 

0,10-0.60 

220 

55,0 

0.436 

0.15-0.65 

OVERSHOOT 

220* 

55,0 

0.44 

0.14-0.90 

NOT  USED 

233 

65.1 

0.309 

0.06-0,65 

234 

80.1 

0.192 

0.06-0.56 

246 

80.0 

0.20 

0.15-0.63 

OVERSHOOT 

245 

89.6 

0.126 

0.10-0.98 

NOT  USED 

247 

89.8 

0.169 

0.14-0.66 

223 

100. 

0  0.151 

0.29-1,00 

OVERSHOOT 

248 

100. 

1  0.135 

0.11-0.64 

OVERSHOOT 

♦HEWLETT-PACKARD  175  OSCILLOSCOPE.  TIME  BASE  CA.  5  ^SEC,  ALL  OTHER  RECORDS  FROM 
TEKTRONIX  454  OSCILLOSCOPE.  TIME  BASE  CA.  l^SEC 
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pressure  at  the  end  of  the  reaction  zone  was  used  as  the  loading 
pressure  of  the  PMMA .  Subsequent  work  has  shown  that  measurements 
at  small  attenuation  are  not  reproducible;  this  is  attributed  to 
the  very  steep  slope  of  the  pressure-time  curve  produced  by  the 
standard  tetryl  donor  of  the  LSGT .  For  this  reason,  all  values 
at  x  <  10  mm  are  considered  nominal  even  in  the  range  x  >.  0.86  mm 
where  measurements  of  u  seem  reproducible.  Moreover,  this  nominal 
range  is  started  (x  =0)  at  a  value  measured  with  a  plane  wave 
boosted  tetryl  loading  of  the  PMMA.  Beyond  x  >.  10  mm,  the  original 
pressure  spike  has  -been  broadened  so  that  its  duration,  seems  suffi¬ 
cient  to  impart  a  particle  velocity  to  the  foil  consistent  with 
the  maximum  amplitude  of  the  pressure  pulse.  In  group  (1), 

Shots  289,  106,  118,  and  120  have  not  been  used. 

In  group  (2),  the  HP175  records  are,  with  one  exception,  from 
the  same  shots  on  which  Tektronix  records  were  made.  Since  we 
consider  the  latter  to  be  in  principle  more  accurate,  the  only  HP 
record  we  have  used  is  that  from  Shot  239  on  which  no  Tektronix 
record  was  obtained.  The  justification  for  using  that  particular 
record  is  the  generally  close  agreement  between  the  HP  and  Tektronix 
results  for  the  same  shot,  as  shown  by  the  data  of  Table  8. 

The  data  selected,  as  above,  from  Table  8  are  plotted  in 
Figure  8a  which  also  shows  two  additional  points  rejected  for  incon¬ 
sistency  with  the  rest  of  the  data.  The  final  calibration  data  have 
been  joined  by  use  of  a  French  curve.  That  was  done  as  a  preliminary 
to  further  study,  and  to  illustrate  the  marked  similarity  between 
the  u-x  curve  of  Figure  8a  and  the  U-x  curve  of  Figure  4  (previous 
section)  in  exhibiting  "humps"  centered  at  x  =  30-35  mm. 

It  is  important  to  note  that  the  data  retained  at  x  =  0  is 
from  Shot  185;  this  used  a  0.001  inch  thick  foil,  a  plane  wave 
booster  and  only  25  mm  of  tetryl.  This  configuration  gives  a  larger 
velocity  (by  0.28  mm/psec)  than  that  using  51  mm  of  tetryl  initiated 
at  a  point;  see  Shot  289  in  Table  8.  The  reason  for  this  is  that 
the  plane  wave  configuration  gives  a  longer  pulse  to  which  the  foil 
can  respond  more  readily.  Thus  the  results  are  more  representative 
of  the  particle  velocity  at  the  interface  than  are  those  from  point 
detonated  charges. 

An  estimate  of  the  precision  of  the  measurements  shown  in 
Figure  8a  can  be  obtained  by  examining  the  replicate  (or  nearly  so) 
shots  listed  in  Table  9,  The  absolute  differences  within  the  sets 
varies  from  0.008  to  0.04  mm/psec.  These  ranges,  as  percentages  of 
the  average  value  of  the  set  of  replicate  shots,  vary  from  0.9  to 
11.2%.  (The  value  of  32%  at  x  *  90  mm  has  been  rejected  as  intol¬ 
erable,  and  the  data  of  Shot  245  have  been  discarded  as  inconsistent 
with  values  observed  at  the  adjacent  stations,  x  «  80  and  x  ®  100  mm. 
See  Figure  8a.)  Because  the  measured  value  u  varies  by  a  factor 
of  10  for  51x1  100  mm,  it  is  not  possible  to  obtain  the  same 
percentage  error  over  the  range.  At  5  mm,  the  shock  is  attenuating 
rapidly  and  the  flow  of  the  detonation  products  frequently  causes 
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Table  9 

RESULTS  FROM  REPLICATE  SHOTS  (TETRYL  DONOR) 


Shot  # 

X 

nun 

Peak  Particle  Velocity 
mm/Msec 

Range 

mm/ysec 

100  x  ran 
average 

241 

5.0 

1.74 

244 

5.0 

1.78 

0.040 

2.3 

258 

5.0 

1.75 

222 

35.0 

0.912 

230 

35.2 

0.904 

0.008 

0.9 

234 

80.1 

0.192 

246 

80.0 

0.200 

0.008 

4.1 

245 

89.6 

0.126* 

247 

89.8 

0.169 

0.043 

32.0* 

223 

100.0 

0.151 

248 

100.1 

0.135 

0.016 

11.2 

♦Rejected  on  basis  of  consistency 
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electrical  noise.  Here,  the  value  of  u  is  relatively  large,  so  that 
a  range  of  0.04  in  the  value  of  u  does  not  make  a  large  percentage 
error.  At  100  mm,  the  shock  is  attenuating  very  little,  and  the 
difference  between  values,  0.016  is  small.  As  a  percentage,  however, 
the  range  is  large  because  of  the  small  magnitude  of  u.  At  this  time, 
we  cannot  measure  values  of  u  as  small  as  0.15  mm/ usee  as  precisely 
as  we  would  like  to. 

Another  approach  to  estimating  the  error  in  an  EMV  gage 
measurement  is  given  in  the  Appendix.  That  estimate  is  for  the 
determination  of  a  single  value  of  u  from  an  oscilloscope  record. 

The  results  are  useful  here,  even  though  wo  have  been  discussing 
peak  values  of  a  which  we  obtain  by  fitting  the  data  to  a  straight 
line  and  extrapolating  back  to  zero  time.  The  estimated  error  ranges 
from  6.2%  for  u  =  0.10  mm/ysec  to  5.6%  for  u  =  2.0  mm/ usee.  The 
slight  variation  of  the  error  with  the  magnitude  of  u  is  due  to  the 
reading  of  the  calibrating  voltage  on  a  three  digit  voltmeter.  With 
a  more  precise  voltmeter,  the  error  would  be  independent  of  the  value 
of  u.  This  comes  about  because  we  change  the  gain  of  the  oscillo¬ 
scope  amplifier  so  that  we  have  about  the  same  deflection  on  the 
record  regardless  of  the  actual  value  of  u  which  is  expected.  This 
is  equivalent  to  changing  scales,  for  example,  on  a  voltmeter,  so 
that  all  readings  are  taken  with  about  the  same  deflection  of  the 
indicator.  This  independence  of  the  estimated  error  of  the  magnitude 
of  the  value  of  u  is  counter  to  our  experience  with  replicate  shots 
as  described  above.  This  disagreement  may  be  due  to  the  fact  that 
our  experimental  arrangements  are  not  reproducible,  as  noted  in  the 
Appendix. 

Figure  8b  shows  the  calibration  data  and  the  data  computed  usina 
WONDY^  for  tetryl  loading  of  PMMA.  Although  the  two  curves  seem 
nearly  the  same  at  small  x  values,  the  experimental  is  definitely 
above  the  computed  curve  in  the  range  15  <.  x  <.  35  mm.  Possible 
reasons  for  this  are  one  or  more  of  the  following: 

a.  Our  extrapolation  procedure  leads  to  too  high  values  of  u. 

b.  The  shock  loading  pressure  (C-J,  not  von  Neumann)  from 
tetryl  was  too  low.* 

c.  The  Hugoniot  used  for  PMMA  was  inaccurate. 

Case  a  alone  is  quite  inadequate  to  explain  the  discrepancy.  The 
velocity  of  the  gage  should  be  equal  to  or  less  than  the  material 
velocity.  Hence  if  we  replace  the  values  obtained  by  extrapolating 
to  zero  time  with  the  maximum  value  of  u  actually  recorded,  it  should 
be  <.  the  true  particle  velocity.  Examination  of  the  records  shows  that 


♦Because  of  zoning  in  the  tetryl,  the  actual  value  of  the  pressure  in 
the  computation  will  be  less  than  the  nominal  C-J  pressure  which  it 
approaches  only  as  the  zone  size  approaches  zero.  In  the  WONDY 
computation  the  zone  size  was  0,05  cm  in  tetryl  and  0.01  cm  in  PMMA. 
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aueh  a  procedure  in  this  range  of  x  lowers  l  he  u  value  by  0,01  mn'^nase 
at  mostj  uhe  difference  wo  wish  to  explain  is  about  tin  this  value. 

Wo  have  already  mentioned  that  04  conditions  were  used  in  the 
computations  instead  of  von  Neumann  conditions,  Hence  the  loading 
shook  pressure  used  (case  b)  wan  too  low,  or  too  much  lowered  by 
the  none  also,  or  both.  The  liugonlot-  data  used  (caae  o)  waa  that 
derived  from  the  previous  1#8ot  calibration,  and  hanee  loan  accurate 
than  those  derived  from  the  :  ©sent  work,  Case  b,«  too  low  an  initial 
pressure,  seems  the  Ukeliesv  single  oaue©  for  the  unexpected  relation, 
u  (experimental)  x  u  (computed), 

We  make  the  assumptions  that  the  Hugoniot,  used  will  give  a  u-x 
curve  of  approximately  the  correct  shape,  and  that  the  shock  loading 
pressure,  as  used  in  the  computations,  was  too  low.  Accordingly, 
the  WONDY  curve  of  Figure  8  b  has  been  raised  (by  0,12*1  mm/i>«e©)  in 
Figure  l)q  so  that  it  approximates  the  experimental  curve  in  the  range 
fl  i  x  £,  15  mm.  This  illustration  strongly  suggests  that  the  on-axis 
flow  is  one-dimensional  for  0  i  x  i  39  mm,  and  becomes  two-dimensional 
with  the  arrival  of  strong  lateral  rarefaction  waves  at;  about 
k  »  33  mm,  This  suggestion  is  reinforced  by  the  data  treatment  which 
will  now  be  described. 

Figure  9a  could  be  used  to  make  graphical  interpolations  between 
data  points,  However,  analytical  relations  which  fit  the  data  arc 
far  more  convenient  for  interpolation,  Moreover,  in  this  case  such 
fits  prove  helpful  in  interpreting  the  data. 

A  polynomial  could  be  used  to  fit  the  data  over  the  entire 
range  of  x,  or  a  set  of  polynomials  could  bo  used  to  fit  subsets 
of  the  data,  (In  this  case,  a  curve  very  like  the  U-x  curve  of 
Figure  4  of  the  previous  section  would  probably  bo  obtained).  Poly¬ 
nomials  have  bad  features,  as  haa  been  related  before.?  In  an  attempt 
to  discover  a  more  appropriate  function,  the  data  were  plotted  on 
semi -log  paper,  On  the  semi -log  plot  shown  in  Figure  9,  the  d  ita 
between  0  and  35  mm  appear  to  be  exponential  in  a  manner  similar 
to  that  found  in  the  radioactive  decay  of  two  isotopes  having  signif¬ 
icantly  different  half-lives.  The  appropriate  function  seems  to  boi 


u  a  A  ©xp(-Bx)  +  C  exp(-Dx)  (8) 


This  function  fits  the  data  very  well)  sea  Table  10,  The  coeffi¬ 
cients,  in  the  order  in  which  they  appear  in  the  above  equation,  and 
their  quadratic  mean  errors  (QME's)  are: 


Coe  f f icient 

m 

1.7342 

0.04  3 

0 .01832 

0.001 

0.6602 

0.041 

0.2794 

0.030 
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Table  10 

FIT  TO  TETRYL  D  V  FOR  x  i  35.2  nun 
FUNCTION  IS  u  =  a  exp(-Bx)+C  exp(-Dx) 


X 

u 

u 

(Calculated) 

Residual 

mm 

mm/psec 

mm/psec 

mm/p sec 

0.00 

2  .  190 

2.39441 

.004 

0.86 

2.250 

2.22600 

-.024 

0.86 

2.220 

2.22600 

.006 

2.00 

2.050 

2.04874 

-.001 

3.00 

1.880 

1.92604 

.046 

5.00 

1.740 

1.74415 

.004 

5.00 

1.780 

1.74415 

-.036 

5.00 

1.750 

1.74415 

-.006 

10.00 

1.490 

1.48141 

-.009 

15.00 

1.300 

1.32354 

-.024 

20.00 

1.190 

1.19983 

-.010 

25.00 

1.110 

1.09206 

-.018 

30.00 

1.010 

0.99506 

-.015 

35.00 

0.892 

0.90694 

.015 

35.20 

0.904 

0.90358 

.000 
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The  QME  of  the  fit  is  0.019.  The  residuals  are  given  in  the  fourth 
column  of  Table  10.  The  largest  residual  is  for  x  =  3.0  mm;  it  is 

0.046,  or  about  2.5%  of  the  measured  value.  A  residual  of  this  size 

is  tolerable  in  view  of  the  range  between  replicates  as  discussed 
above,  and  also  that  x  =  3  ram  is  within  the  range  of  "nominal"  cali¬ 
bration.  Outside  this  range  of  u,  (i.e.,  x  <  10  mm)  the  largest 
residual,  on  a  percentage  basis,  is  1.8%. 

There  appears  to  be  a  distinct  change  in  the  slope  of  the  u  vs  x 
relation  at  about  35  mm,  see  Figure  8c.  Data  beyond  35  mm  must  be 
fitted  with  another  function.  As  shown  in  the  semi-log  plot  of 
Figure  9,  two  straight  lines  can  be  drawn  through  the  data;  the  first 

serves  for  35  <_  x  <_  75  mm,  while  the  second  serves  for  75  £.  x  <.  100. 

Thus  it  appears,  at  first  glance,  that  two  or  more  functions  are 
required,  and  that  a  distinct  change  in  slope  might  be  introduced 
at  75  mm.  An  examination  of  Figure  8a  shows  that  a  discontinuity  in 
the  slope  of  the  u-x  curve  at  75  mm  is  not  evident.  Moreover,  it  is 
not  necessary  to  separate  these  data  into  two  subsets  for  x  >  35  mm. 

The  strategy  is  to  fit  the  quantity  (u  -  A)  vs  x.  That  is,  subtracting 
a  judiciously  chosen  value  of  A  from  the  values  of  u  yields  a  straight 
line  plot  on  semi-log  paper.  The  results  of  such  a  fit  are  given  in 
Table  11.  The  coefficients  and  their  QME's  for  the  relation 


u  =  A  +  B  exp(-Cx)  (9) 

are 

Coefficient  QME 

0.0921  0.0065 

3.7038  0.135 

0.0435  0.0011 

The  QME  of  the  fit,  0,008,  is  about  half  of  the  QME  of  the  fit 
given  in  Table  10.  Examination  of  the  residuals  shows  that  the  fit 
is  really  very  good,  considering  that  the  function  has  only  3  param¬ 
eters.  Unavoidably,  the  residuals  for  very  low  u  values  are  a  high 
percentage  of  u. 

It  must  be  noted  in  passing  that  some  of  the  data  shown  in 
Table  11  are  from  shots  using  pentolite  donors.  In  a  later  section, 
it  will  be  noted  that  tetryl  and  pentolite  give  values  of  particle 
velocity  which  are  not  distinguishable  at  high  attenuations.  Hence, 
the  u  vs  x  curves  are  allowed  to  converge  for  some  value  of  x  less 
than  80  mm. 

Note  that  in  Table  11,  two  or  three  values  of  u  are  given  for 
several  individual  values  of  x.  If  all  values  of  u  were  given  the 
same  weight,  the  function  would  have  fitted  these  multivalued  points 
more  closely  than  the  single  value  points.  This  was  prevented  by 
giving  the  multivalued  points  a  weight  of  1.0,  while  the  single 
valued  points  were  given  a  weight  of  2.0. 
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Table  11 


FIT 

TO  TETRYL  DATA 

FOR  x  i.  35  mm 

AND 

PENTOLITE  DATA 

FOR  x  >.  80  mm 

THE 

FUNCTION  IS  U  = 

'  A+B  exp(-Cx) 

X 

u 

Weight 

u 

Residual 

mm 

mm/ysec 

(Calculated) 

mm/usec 

mm/ysec 

35.0 

. 89200 

1.0 

.90010 

.008 

35.2 

.90400 

1.0 

.89310 

-.011 

39.9 

.73500 

2.0 

.74499 

.010 

45.0 

.62500 

2.0 

.61509 

-.010 

55.0 

.43600 

2.0 

.43061 

-.005 

65.1 

.30900 

2.0 

.31025 

.001 

80.0 

.20000 

1.0 

.20620 

.006 

80.0 

.19400 

1.0 

.20620 

.012 

80.1 

.19200 

1.0 

.20571 

.014 

89.8 

.16900 

1.0 

.16660 

-.002 

100.0 

.15100 

1.0 

.13991 

-.011 

100.0 

.14100 

1.0 

.13991 

-.001 

100.1 

. 13500 

1.0 

.13970 

.005 

101.6 

.14800 

1.0 

.13670 

-.011 
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Figure  10  shows  the  u-x  data  plotted  along  with  the  solid  lines 
of  the  fitted  equations  (Equations  (8)  and  (9)").  This  treatment 
has  resolved  the  "hump"  of  Figure  Ra  into  a  cusp  between  the  two 
analytic  curves  near  x  =  35  mm,  at  about  the  location  of  the  maximum 
amplitude  of  the  so-called  hump.  Figure  10  illustrates  much  more 
clearly  than  Figure  8c,  the  termination  of  one-dimensional  flow  by 
the  arrival  of  a  lateral  rarefaction  at  the  axis. 

Not  enough  measurements  were  made  in  the  earlier  work  to 
reveal  the  cusp  which  is  evident  in  the  new  data  at  about  35  mm. 

It  appears  that  the  slope  of  the  u  vs  x  curve  changes  abruptly  at 
this  distance.  This  is  doubtlessly  caused  by  the  interaction  of 
rarefaction  waves  which  originate  at  the  periphery  of  the  PUMA 
cylinder  at  the  interface  between  it  and  the  explosive.  The  strength 
of  this  rarefaction  increases  because  of  radial  convergence.  When 
it  reaches  the  axis  of  the  PMMA  cylinder,  it  causes  an  abrupt  change 
in  the  flow  due  not  only  to  its  strength  but  also  to  the  gradient 
of  flow  variables  across  it.  Another  manifestation  of  this  wave 
interaction  was  pointed  out  in  a  previous  section  where  the  shock 
velocity  was  measured  as  a  function  of  x.  Indeed,  the  similarities 
of  Figures  4  (previous  section)  and  8a  indicate  that  a  similar 
analytical  treatment  of  the  U-x  data  might  result  in  showing  a  cusp 
at  x  =  35  mm  in  the  shock  velocity  data  also. 

Kamegai  and  Erkman-^  used  a  two-dimensional  elastic-elastic  code, 
HEMP,  to  study  the  diverging  shock  wave  in  the  PMMA  in  the  LSGT 
configuration.  This  earlier  work  was  done  without  knowledge  of  the 
particle  velocity  measurements  discussed  above.  Two  different  grid 
sizes,  0.06  cm  sq.  and  0.12  cm  sq.,  were  used  in  the  2-D  work. 
Unfortunately,  only  the  0.12  cm  sq .  zone  results  covers  the  range 
0  to  100  mm.  As  shown  in  Figure  11,  the  2-D  study  reproduces  the 
shaoe  of  the  u  vs  x  curve  obtained  by  the  EMV  method  but  the  code 
results  are  larger.  These  larger  values  of  u  may  bo  due  to  the  use 
of  a  smaller  value  of  y  (and  correspondingly  higher  value  of  P<gj) 
for  the  explosive  gases?  a  Y  of  2.54  was  used  while  experimental 
work^  yields  a  y  of  3.00.  With  this  in  mind,  the  calculated  and 
experimental  u  vs  x  curves  agree  well.  Shifting  the  computed  curve 
down  about  0.07  mm/psec  gives  better  agreement,  of  course,  and,  more 
important,  clearly  illustrates  the  same  type  of  abrupt  change  in 
slope  at  ■-  35  mm  in  the  computed  curve.  An  analytical,  treatment 
similar  to  that  used  on  the  experimental  data  might  well  demonstrate 
a  cusp  in  the  2-D  curve,*  and  appears  to  be  a  more  reliable  way  to 
locate  the  arrival  of  the  lateral  rarefaction  than  the  method  chosen 
in  reference  (10).  That  method  was  based  on  finding  where  the 
calculated  shock  front  deviated  from  spherical  expansion.  The  result 
was  that  the  flow  seemed  to  be  spherical  out  to  s  50  mm  whore  the 
side  rarefaction  finally  reached  the  axis  of  the  cylinder.  It  now 


*Of  course,  2-D  computations  of  more  closely  spaced  points  should 
do  the  same  thing. 
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appears  that  the  front  is  affected  by  lateral  rarefactions  at  about 
35  mm;  beyond  that  distance  the  shock  is  two-dimensional. 

In  previous  work,  velocities  of  the  free  surfaces  of  the  ends 
of  cylinders  such  as  shown  in  Figure  5a  were  measured.  If  the 
shocked  PMMA  behaves  as  an  inviscid  fluid,  the  particle  velocity 
is  half  the  measured  free  surface  velocity.  Figure  12  shows  the 
difference  between  this  estimate  and  the  measured  particle  velocities 
in  PMMA.  Representative  values  are  as  follows: 


Gap  Length 
mm 


_  u,  mm/tjsec _ 

Free  Surface  Method  EMV  Method 


Difference 
mm/ y sec 


1.09 

0.54 

0.21 

0.145 


1.19 

0.63 

0.22 

0.151 


0.10 

0.09 

0.01 

0.006 


The  free  surface  method  gives  significantly  smaller  values  of  u  at 
20  and  45  mm.  At  75  and  100  mm  the  two  methods  are  in  substantial 
agreement.  Thus  the  new  work  shows  significantly  larger  particle 
velocities  of  moderate  (and  presumably,  low)  attenuations;  this  will 
result  in  larger  pressure  also.  It  is  chiefly  the  improvement  in 
values  of  u  that  requires  the  present  recalibration  of  both  the 
tetryl  and  the  pentolite  loaded  PMMA. 

Data  from  shots  using  pentolite  donors  are  given  in  Table  12 
and  in  Figure  13.  Table  12  has  the  same  format  as  Table  8  for  ease 
in  comparing  results.  Another  aid  in  comparing  the  two  sets  of 
data  is  the  curve  in  Figure  13  which  represents  the  fits  to  the 
tetry.l  data.  For  10  <l  x  <.  65  mm,  pentolite  donors  give  slightly 
greater  values  of  u.  This  is  to  be  expected  because  pentolite  has 
an  estimated  Chapman-Jouguet  pressure  of  216  kbar  while  tetryl  has 
a  CJ  pressure  of  190  kbar.H  This  difference,  for  some  unknown 
reason,  is  not  apparent  in  our  data  for  x  <  10  mm. 

It  should  be  remarked  at  this  point  that  for  any  x  >  10  ram, 
the  difference  between  u  measured  for  pentolite  and  tetryl  loading 
is  experimentally  insignificant.  This  is  indicated  both  by  our 
error  analysis  and  by  the  range  found  for  replicate  measurements. 
However,  most  pentolite  data  were  consistently  above  the  analogous 
tetryl  data.  Hence  we  have  presented  two  distinct  calibration  (u-x) 
curves  with  the  small  differences  between  them  shown  in  Figure  13. 

Again,  we  report  results  from  the  Hewlitt  Packard  oscilloscope 
but  do  not  use  them  in  the  subsequent  curve  fitting.  Results  of 
Shot  260  were  not  used  because  of  the  use  of  a  "radial  foil"  (i.e., 
the  leads  from  the  pickup  came  out  of  the  PMMA  cylinder  parallel  to 
a  radius  rather  than  as  shown  in  Figure  5a) .  Shot  261  was  rejected 
because  it  gave  results  that  were  less  than  those  in  tetryl  shocked 
PMMA  at  the  same  position.  The  same  is  true  of  Shot  210. 
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TABLE  12 


H 


1 


1 


i 


PEAK  PARTICLE  VELOCITIES  FOR  PENTOLITE  DONOR 


COMMENTS 


SHOT 

X 

U 

TIME  INTERV 

NO 

MM 

MM//XSEC 

FOR  LS  FIT 

286 

0.0 

2.37 

0.06-0.17 

277 

5.0 

1.80 

0,03-0,45 

278 

5.0 

1.74 

0.018-0.96 

275 

10.0 

1.494 

0.05-0.54 

261 

10.0 

1.39 

0.12-0.63 

262 

10.0 

1.51 

0.09-0.35 

276 

15.0 

1.351 

0.10-0.35 

263 

20.0 

1.257 

0.13-0.63 

264 

20.0 

1.22 

0.04-0.39 

270 

25.0 

1.12 

0.13-0,65 

260 

25.0 

1.09 

0.10-0.60 

271 

25.0 

1.11 

0.05-0.56 

235 

30.0 

1.038 

0.06-0.45 

2  35* 

30.0 

1.05 

0.20-1.50 

265 

30.0 

1.01 

0.07-0.56 

225 

35.0 

0.942 

0.05-0.55 

225* 

35.0 

0.945 

0,10-0.85 

229 

38.1 

0.826 

0.05-0.55 

236 

?9.9 

0.786 

0.03-0.47 

236* 

39.9 

0.78  7 

0.20-1.00 

219 

45.0 

0.661 

0.09-0.60 

210 

45.0 

0.613 

0.10-1,10 

211 

54.9 

0.431 

0.27-1.09 

221 

55.0 

0.451 

0,07-0,60 

237 

65.1 

0.313 

0,03-1.00 

266 

80.0 

0.  194 

0, 09-0.54 

224 

100.0 

0,141 

0.09-0.60 

194 

101.6 

0.149 

0. 10-0. 60 

NOT  USED*  LOWER  THAN  TETRYL 


OVERSHOOT 

RADIAL  FOIL*  NOT  USED 
TYPICAL  RECORD 
TYPICAL  RECORD 
NOT  USED 


NOT  USED 

TEST  WITH  TWO  PIECES  OF  PMMA  •  NOT  USED 
NOT  USED 

LESS  THAN  TETRYL  AT  SAME  DISTANCE 
EYE  FIT  IS  USED 


^  °™"  "EC°R‘’S  FR°" 


PEAK  PARTICLE  VELOCITY,  mm 
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FIG,  13  EXPERIMENTAL  RESULTS  (PENTOLITE  DONOR)  COMPARED  WITH  RESULTS 
USING  TETRYL  DONOR 
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Our  experience  with  t:ho  ealibrat ton  for  tetryl  showed  that  w® 
could  expect  erratic  results  close  to,  and  at,  the  interface  %  At 
the  interface,  using  a  51  mm  long  point,  initiated  ptmtolite  uharye, 

we  measured  a  particle  velocity  of  2,37  mm/imuo,  We  axp«ot  that, 

results  from  a  plane  wave  boosted  (PWB)  charge,  uainq  25  mm  of  pente- 

l.ite,  would  be  greater  than  2.37,  as  it  was  for  totrylj  tui®  above. 

The  difference  for  tetryl  was  0.28  mm/usoo.  Wt>  add  this  increment 
to  the  value  from  Shot  286,  Tabic  12,  to  obtain  2.65  mm/ usee  as  our 
interface  value.  As  will  be  shown  later,  this  gives  u»  a  reasonable 
calibration  in  the  nominal  region,  0  ^  x  10  mm. 

The  pentolite  data  are  fitted  with  the  double  exponential, 
Equation  (8),  over  the  first:  35  nun.  The  coefficients  and  their  QME's 
are 


Coefficient  QMK 


1,773  0.032 

0.0184  0.0008 

0.876  0.036 

0.349  0.032 


Residuals,  which  are  given  in  Table  13,  are  all  leas  than  0.03  in 
absolute  value.  The  results  of  the  fit  are  displayed  in  Figure  13 
along  with  the  fit  to  the  tetryl  data,  and  the  experimental  pentolite 
data.  The  fit  appears  to  be  adequate  for  interpolating  over  this 
part  of  the  experimental  data. 

Data  for  pentolite  at  x  '•  35  mm  are  fitted  with  Equation  (9) 
for  which  the  coefficients,  in  order,  and  their  QMK *  s  are 


Coefficient 


QMK 


0.0905  0.0075 
4.0877  0.1588 
0.0445  0,0012 


The  data  and  residuals  are  given  in  Table  14.  In  obtaining  this  fit, 
data  for  tetryl  shots  were  included  for  x  >.  80.0  mm.  Results  of  this 
fit  are  also  shown  i.n  Figure  13.  Note  that  the  results  for  tetryl 
and  pentolite  converge  for  x  >  65  nun.  The  slight  differences  in  the 
two  curves  for  x  >  65  mm  has  no  practical  significance. 

The  two  curves  which  represent  the  tetryl  data  cross  at 
x  =  34.65  mm;  those  representing  pentolite,  at  36  mm.  When  the  two 
u,x  curves  are  plotted  on  a  very  large  scale  and  compared  in  the 
region  33  i  x  i  38  mm,  they  are  practically  coincident  after  one  is 
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Tall  j  0  n 

FIT  TO  1'KNTOl.I  TK  DATA  X  4,  36 ,  0  nun 

TUN  FUNCTION  ffl  u*A  *xp (-&K* *C  axp(-Gx) 
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shifted  1,35  nun,  the  difference  found  in  location  of  the  rarefaction 
waves  for  the  two  boosters,  The  fact  that  these  locations  differ 
has  an  interesting  effect  in  our  pressure  calibration)  this  will  be 
described  in  a  later  section, 

in  summary,  the  peak  particle  velocity  in  the  PMMA  attenuator 
of  the  UJGT  for  tetryl  donor  charges  is  given  by 


u  »  1,7342  exp(-Q. 01852  x)  +  0,6602  axp(-0.2794  x) 


for  x  i,  34, 6 S  mm 
and 

u  *»  0.0021  +  3.7038  exp(-Q.Q435  x)  for  x  >  34. GS  mm. 

For  pentolite  donor  charges,  the  results  are 

u  ®  1.7735  exp(-0. 01Q41  x)  ♦  0.8765  exp(-0.3495  x) 
for  x  <.  36  mm 

and 

u  «>  0,0905  +  4.0877  exp  (-0 .04451  x)  for  x  >  36  mm. 


(10) 

(ID 


(12) 

(13) 


l 


5.  HUGONXOT  OF  PMMA 

The  end  result  of  a  gap  calibration  is  a  curve  relating  pressure 
to  length  of  gap.  The  data  presented  in  the  two  preceding  sections 
could  bo  used  to  provide  the  calibration.  That  is,  for  every  value 
of  x,  we  have  values  of  u  and  U  which  define  a  value  of  P.  Our  data, 
however,  are  not  entirely  accurate  as  has  been  noted  in  previous 
sections.  One  way  to  examine  the  accuracy  of  the  data  is  to  study 
the  Hugoniot  curve  which  results  by  plotting  U  as  a  function  of  u. 

Then  our  results  can  be  compared  with  results  previously  reported. 

In  the  following  discussion,  we  will  be  interested  in  pressures 
up  to  about  200  kbar.  This  is  well  below  the  phase  change  in  PMMA 
which  takes  place  between  212  and  276  kbar. 5 

Data  on  the  Hugoniot  of  PMMA  are  available  in  the  LLL  compendium. 12 
Most  of  the  data  that  were  available  in  1965  were  reviewed  by  Deal 
who  noted  that  there  was  a  large  amount  of  scatter  in  the  data.*3  He 
questioned  the  reproducibility  of  the  material  used  in  the  experi¬ 
ments  —  few  authors  gave  a  thorough  characterization  of  the  plastic 
used.  Some  of  the  scatter  may  well  be  due  to  differences  in  material. 
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Different  techniques  may  also  have  contributed  to  the  scatter.  Since 
1965,  data  have  been  acquired  with  now  techniques  on  easily  identified 
materials.  The  most  aor urate  of  those  uses  optical  interferometry. 

Data  taken  wi^h  thin  technique  at  Sandia  Laboratories  in  the  range  of 
1-60  kb»r^,l.>  are  used  in  the  following  as  a  check  on  our  data. 

Results  from  General  Motors  confirm  the  low  pressure  results  reported 
by  Sandia,  Both  sots  of  data  are  from  samples  cut  from  sheet  stock 
of  Plexiglas,  produced  by  Rohm  and  Haas  and  designated  Plexi¬ 
glas  II  UVA,  an  ultraviolet  absorbing  grade.  This  is  the  same  material 
which  we  use,  the  exception  being  that  we  use  bar  stock.  As  noted 
above,  optical  interferometry  was  usod  to  study  shocks  in  PMMA  at 
Sandia,  At  General  Motors,  the  data  was  acquired  by  the  use  of 
transducers,  a  method  which  .is  fairly  new  in  this  type  of  work. 

Because  the  data  from  Sandia  extends  over  a  greater  range  of  pressures, 
we  depend  on  them  in  the  following  discussion. 

In  these  recent,  works,  it  has  been  shown  that  PMMA  behaves  in  a 
complicated  way  when  stressed  by  shock  waves.  It  is  generally  conceded 
that  the  material  is  viscoelastic,  and  exhibits  strain  rate  effects „ 
There  is  also  evidence  that  it  behaves  elastoplastically ,  at  least 
when  shocked  to  states  above  25  kbar.  One  observation  in  refer¬ 
ence  (14)  is  that  tho  shock  .velocity  in  an  impaetod  piece  of  PMMA 
depends  on  the  thickness  of  the  piece.  That  is,  even  though  there 
is  no  interaction  of  waves  from  the  back  of  the  projectile  with  the 
shock  front  in  tho  sample,  tho  shock  velocity  changes  with  distance. 
This  is  consistent  with  the  observation  that  PMMA  is  a  viscoelastic 
solid.  Another  observation  is  that,  for  example,  impact  by  a  PMMA 
projectile  at  a  velocity  of  about  0.6  mm/nsoc  produces  an  abrupt  rise 
of  particle  velocity  to  about  0.20  mm/usec,  followed  by  a  slower  rise 
to  about  0.30  mm/usec.  The  latter  is  equivalent  to  a  pressure  of 
about  12.3  kbar.  This  same  typo  of  behavior  is  observed  at  19.4  kbar, 
while  at  60  kbar  the  entire  loading  event  is  abrupt  as  it  would  be  in 
a  metal. 15  At  low  pressures,  we  really  should  not  use  the  words 
"shock"  and  "Hugoniot"  when  we  describe  the  response  of  PMMA  when  it 
is  impacted  by  a  projectile,  or  loaded  by  an  explosive.  For  conven¬ 
ience,  however,  we  will  continue  to  use  the  familiar  nomenclature. 

When  we  use  the  term  "shock  wave",  we  mean  the  entire  loading  process 
which  takes  place  in  PMMA  whether  or  not  it  is  entirely  abrupt. 

Our  experimental  relation  between  U  and  u  is  shown  in  Figure  14. 
This  is  a  graph  of  the  data  presented  in  the  previous  two  sections  of 
this  report  for  pentolite  donor  charges.  Each  point  in  the  graph 
corresponds  to  a  value  of  x,  the  distance  from  the  donor  explosive. 

Some  representative  values  of  x  are  indicated  on  the  graph  as  an 
aid  in  the  following  discussion.  Our  results  are  not  linear,  even 
at  the  high  pressure  end  of  our  set  of  data.  A  consequence  of  this 
is  that  extrapolation  is  more  of  a  problem  than  it  would  be  if  our 
data  were  linear  at  the  higher  pressures.  A  procedure  for  extrap¬ 
olating  will  be  discussed  below.  Curvature  at  low  pressure  is 
expected.  Data  acquired  previously  at  NOL  showed  that  PMMA  had  an 
anomolous  Hugoniot. 4  Data  from  SRI,^7  Sandia  Laboratories ,1 5  ancj 
General  Motors^  confirm  the  observation  of  this  behavior. 
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Schuler  and  Nunziato*-4  have  an  expression  which  fits  their  data 
and  that  of  Barker  and  Hollenbach. *3  They  chose  to  relate  the  stress 
to  the  strain  rather  than  to  relate  shock  and  particle  velocities. 
Their  expression  is 


c  *  87 . 066e  +  858. 71e2  -  7063.23 e3  +  22040. 8e4  (14) 


where  a  is  the  stress  in  kbars  and  e  is  the  strain  (engineering) . 
Using  the  relation 


e  =  1  -  VA0  (15) 

where  V  is  the  specific  volume  and  Vg  is  the  specific  volume  at 
ambient  stress,  and  the  Hugoniot  relations,  there  results, 


e  =  u/U  (16) 

o  =  p0uU  (17) 


The  above  relations  permit  the  construction  of  a  u,U  curve  from  the 
Sandia  relation,  Equation  (14).  With  this  representation,  we  have  a 
continuous  curve  covering  about  60  kbar  of  the  PMMA  Hugoniot.  Curva¬ 
ture  at  low  pressure  is  more  pronounced  than  it  is  in  our  data,  see 
Figure  14.  This  curvature  is  substantiated  by  the  results  from 
General  Motors^  which  are  also  shown  in  Figure  14.  Schuler  and 
Nunziatol^  did  not  use  the  two  points  from  Barker  and  Hollenbachi4 
for  which  u  is  about  0.61  and  0.625  mm/usec  in  their  determination  of 
the  constants  in  Equation  (14).  Their  reason  for  not  usinq  these  data 
is  that  they  were  obtained  from  experiments  in  which  the  projectiles 
were  tungsten  carbide.  These  experiments  could  not  be  analysed  as 
satisfactorily  as  the  other  experiments  in  which  the  projectiles  were 
PMMA.  Hence  they  decided  to  ignore  these  two  points  rather  than  to 
take  a  chance  on  less  than  satisfactory  results. 

The  curve  shown  in  Figure  14  for  the  Schuler  and  Nunziatol5  data 
was  constructed  from  "equilibrium"  values  of  particle  velocity.  In 
some  instances,  a  time  interval  as  long  as  one  usee  was  required  for 
the  instantaneous  u  to  build  up  to  its  higher  equilibrium  value. 

Since  this  build-up  time  is  often  greater  than  our  time  of  measurement 
of  u,  there  is  some  question  of  the  applicability  of  the  equilibrium 
curve  to  our  data.  To  be  sure,  when  the  Liddiard4  data  (U,  1/2  UfS) 
are  corrected  to  U,  u  data  by  means  of  Figure  12,  they  fall  on  the 
equilibrium  curve.  However,  this  may  be  a  fortuitous  agreement. 

The  difference  between  the  curve  derived  from  Equation  (14)  and 
the  curve  which  represents  our  data  is  small.  This  is  shown  in 
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Figure  15  where  the  difference  in  pressure  between  the  two  represen¬ 
tations  is  plotted  as  a  function  of  x.  In  constructing  this  plot,  we 
used  the  U  vs  x  data  reported  in  Section  3  and  the  u  vs  x  of  Section  4 
of  this  report  for  pentolite  donors.  Pressure,  based  on  the  Sandia 
curve,  is  also  shown  on  the  plot.  The  difference  in  pressure  is  less 
than  a  kbar  except  for  x  <  13  mm.  Hence  either  representation  could 
be  used  in  our  calibration  over  the  range  of  pressure  up  to  about 
80  kbar. 

We  need  to  extend  our  Hugoniot  representation  to  about  200  kbar 
so  we  can  obtain  a  nominal  calibration  for  the  LSGT  for  small  gaps. 
Examining  Figure  14  shows  that  we  could  use  the  Sandia  curve  for  the 
extrapolation  —  it  reproduces  our  data  to  within  experimental  error 
up  to  a  particle  velocity  of  1.5  mm/ysec  (90  kbar).  This  agreement 
may  be  fortuitious  as  will  be  demonstrated  in  the  following.  Some  of 
the  data  reported  by  Coleburnl^  are  shown  in  Figure  14 ;  these  indicate 
that  the  extrapolated  Sandia  curve  is  too  steep.  The  slope  of  the 
curve  is  about  1.9  for  0.9  <.  u  <.  1.9  mm/ysec.  Coleburn  fitted  his 
data  with  a  quadratic  which  has  a  slope  of  about  1.47  at  u  =  1.5. 
Values  of  the  slope  of  straight  line  fits  to  PMMA  Hugoniot  data  in 
this  range  have  been  reported  as  1.61^'^  and  1.69.4  It  was  decided 
not  to  risk  using  the  extrapolation  with  the  steep  slope,  so  we  tried 
finding  a  straight  line  through  a  part  of  our  data  which  could  be  used 
for  extrapolation.  Placing  a  straight  edge  on  the  plot  shows  that  a 
straight  line  can  be  drawn  through  two  subsets  of  the  data.  These 
subsets  are  for  0.595  <.  u  <.  0.749  and  1.002  <.  u  <.  1.324  .  Ignoring 
the  non-linearity  of  the  intermediate  data  is  justified  as  follows. 
Those  data  which  sag  below  the  straight  line  come  from  measurements 
for  x  between  31  and  41  mm.  It  must  be  recalled  that  the  u  vs  x  curve 
was  treated  so  that  a  cusp  appeared  at  about  35  mm.  If  wave  inter¬ 
action  in  the  cylinder  of  PMMA  can  produce  this  anomolous  behavior  in 
the  particle  velocity  curve,  it  can  surely  produce  similar  behavior 
in  the  shock  velocity.  We  may  have  smoothed  the  x,t  data  and  the 
U,x  data  so  that  we  suppressed  some  of  the  effects  of  wave  interaction 
on  the  results.  We  do  observe  a  hump  in  the  U,x  curve  at  about  35  mm. 
Treated  differently,  this  hump  might  look  more  like  the  cusp  which 
appears  in  our  treatment  of  the  particle  velocity.  The  results  might 
also  eliminate  the  sag  in  our  U,u  curve.  Furthermore,  the  greatest 
divergences  in  shock  velocities  were  observed  at  35  and  40  mm,  see 
Table  6.  These  considerations  do  not  prove  that  those  points  which 
we  ignore  are  really  in  error  —  they  do  point  to  the  possibility. 

The  data  at  the  top  end  of  our  curve  (u  >  1.324)  come  from 
locations  within  16  mm  of  the  donor  explosive.  In  this  region  the 
shock  velocity  changes  rapidly,  with  the  result  that  the  streak  camera 
records  are  difficult  to  differentiate.  One  of  the  two  records 
gave  very  poor  data  over  the  first  8  mm.  Here  too,  a  reasonable  case 
is  available  for  ignoring  the  nonlinear  data  at  ths  high  pressure  end 
of  our  curve.  The  fact  that  nonlinearity  in  this  region  of  the  U,u 
curve  has  not  been  reported  by  other  investigators  has  also  been 
considered. 
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The  line  mentioned  above  is  represented  by 


U  =  2.561  +  1.595  u. 


(18) 


In  previous  work  with  tetryl  donor  charges  the  coefficients  were 
found  to  be  2.57  and  1.61  respectively . 3  Considering  all  the  possi¬ 
bilities  of  errors  in  measurements  and  in  interpretation,  the  two 
results  are  in  good  agreement. 

At  low  pressures  (low  u  and  U) ,  we  are  at  the  other  extreme  of 
our  experimental  measurements.  There  seems  little  doubt  that  the 
Sandial4,15  measurements  should  be  more  accurate  than  any  we  can  make. 
Consequently  we  have  fitted  the  Sandia  data,  shown  in  Table  15,  with 
a  cubic  constrained  to  become  tangent  to  the  straight  line  of  Equa¬ 
tion  (18).  Thus  our  final  selection  for  the  Hugoniot  is 

U  =  2.7228  +  4.0667u  -  10.9051u2  +  10.6912u3,  0.03  £.  u  <  0.5363  (19) 

and 


U  =  2.561  +  1.595u,  u  >  0.5363. 


(18') 


Figure  16  shows  Equations  (19)  and  (18')  as  the  solid  line  curve, 
the  Sandia  data  from  which  Equation  (19)  was  obtained,  and  the  low 
pressure  U,u  data  of  the  present  work.  The  latter,  discarded  in 
selecting  the  Hugoniot,  lead  to  differences  of  <  0.3  kbar  (at  the 
same  x  or  u  value)  in  the  pressure  computed  with  the  Hugoniot  of 
Equation  (19)  and  those  computed  with  the  discarded  U,u  data  (see 
Figure  15) . 


6.  PRESSURE-DISTANCE  CALIBRATION  FOR  THE  LSGT 

We  now  have  the  particle  velocity  vs  distance  relation  for  the 
PMMA  attenuator  of  the  LSGT  in  analytical  form  for  both  tetryl  and 
pentolite  donor  charges  (Section  4) .  These  constitute  our  primary 
calibration  curves  for  reasons  presented  in  the  error  analysis.  In 
addition,  we  have  analytical  expressions  relating  the  particle  and 
shock  velocities  along  the  PMMA  Hugoniot  (Equations  (18')  and  (19)  of 
the  previous  section) .  From  the  above  combination  we  can  calculate 
the  pressure  in  the  attenuator  at  any  distance  from  either  donor 
charge  by  using 


P  =  pg  u  U  (20) 


where  pq  is  the  density  of  PMMA,  taken  to  be  1.185  g/cc.  Table  16 
gives  the  computed  pressures  at  5-10  mm  intervals  in  x. 
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Table  15 

CUBIC  FIT  TO  LOW  PRESSURE 
HUGONIOT  DATA  FOR  PMMA 
(SANDIA  CORPORATION  DATA) 


U 


u 

U 

(Calculated) 

Residual 

mm/ y sec 

mm/ y sec 

mm/ysec 

mm/ysec 

.03045 

2.843 

2.83685 

-0.006 

.03066 

2.844 

2.83757 

-0.006 

.07555 

2.968 

2.97243 

0.004 

.07580 

2.959 

2.97308 

0.014 

.1542 

3.127 

3.12981 

0.003 

.1545 

3.130 

3.13025 

0.000 

.2250 

3.199 

3.20754 

0.009 

.3206 

3.268 

3.25804 

-0.009 

.3196 

3.268 

3.25767 

-0.010 

.4801 

3.349 

3.34476 

-0.004 

.4805 

3.342 

3.34516 

0.003 

.228 

3.22 

3.20986 

-0.010 

.500 

3.32 

3.36629 

-0.046 

NOTE:  Last 

two  points 

from  Schuler  and  Nunziato; 

remainder 

from  Barker  and  Hollenbach. 
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Since  the  unit  of  length  used  in  the  gap  test  is  the  thickness 
of  a  card  (0,01  inch),  it  is  convenient  to  tabulate  the  pressure  at 
one  card  intervals.  This  can  easily  be  done  by  use  of  a  digital 
computer;  Tables  17  and  18  contain  the  final  calibration  values  for 
tetryl  and  pentolite  donors,  respectively.  The  gap  length  in  these 
tables  is  the  sum  of  the  number  in  the  left  column  and  the  number 
at  the  top  of  a  column.  Each  table  gives  P  for  zero  through  400 
cards.  Although  pressure  is  reported  to  0.1  kbar,  the  fractional 
part  is  not  significant. 

A  rigorous  analysis  of  the  possible  error  in  the  pressure  has  not 
been  carried  out.  An  analysis  of  the  error  expected  in  the  measured 
particle  velocity  gives  an  estimate  of  ±6%  in  an  individual  determi¬ 
nation  (see  Appendix) .  Replicate  measurements  over  the  range  of 
0.2-1. 5  mm/ysec  show  differences  of  4%  or  less  and,  therefore,  suggest 
the  ±6%  is  probably  pessimistic;  however,  replicates  at  0.15  mm/ysec 
show  a  difference  of  11%.  The  error  in  P  must  be  at  least  as  large 
as  the  error  in  u;  it  is  probably  somewhat  larger  because  of  any  error 
in  the  Hugoniot  we  have  selected. 

As  in  the  case  of  particle  velocity,  at  any  given  x  >.  10  mm 
(39  cards),  the  difference  in  pressure,  obtained  by  tetryl  and 
pentolite  loading,  is  not  experimentally  significant.  The  numerical 
difference  (values  of  Table  18  minus  those  of  Table  17)  runs  from 
32  kbar  at  x  =  0  cards  to  0  at  400  cards.  The  maximum  percentage 
difference,  in  the  range  x  >  10  cards,  is  about  7%  (2 . 6"  kbar)  ancl 
occurs  at  x  of  140-150  cards.  This  is  the  direct  result  of  a  similar 
maximum  difference  in  particle  velocity  caused  by  differences  in  the 
location  of  the  rarefaction  waves  (cusps  in  the  fitted  u-x  curves)  in 
the  two  cases.*  It  was  remarked  in  Section  4,  that  shifting  one  of 
the  u-x  curves  by  1.35  mm  (or  about  5  cards)  resulted  in  coincidence 
of  the  two  curves.  In  agreement  with  this.  Tables  17  and  18  show 
that  the  pressure  at  x  cards  for  tetryl  is  the  same  (±0.1  kbar)  as 
that  at  (x  +  5)  cards  for  pentolite  in  the  range  of  70-230  cards.  At 
greater  attenuations,  of  course,  the  pressure  is  the  same  (±0.1  kbar) 
for  the  same  number  of  cards;  this  is  true  for  the  range  260-400  cards. 

In  earlier  work,^  two  pairs  of  tests  were  carried  out  in  the 
region  70-230  cards.  For  the  same  acceptor  and  different  donors 
(tetryl  and  pentolite),  the  pairs  showed  differences  of  one  and  zero 
cards  in  the  two  respective  50%  gaps.  This  is  not  the  5  card  differ¬ 
ence  expected,  but  it  indicates  the  same  pressure  amplitude  within 
the  expected  error.  The  expected  error  is  just  too  large  to  resolve 
the  effect  of  different  pressure  profiles  of  the  same  amplitude. 


*Use  of  our  U,u  data  or  the  Sandia  equilibrium  U,u  data  instead  of 
the  Hugoniot  we  have  chosen  has  no  effect  on  the  size  of  this  differ¬ 
ence.  The  .Oandia  data  result  in  lower  absolute  pressure  values;  at 
x  2  80  cards,  the  maximum  difference  is  1.5  kbar  and  1  kbar, 
respectively,  for  tetryl  and  pentolite. 
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Table  1*)  haa  buun  prepared  ««  that  the  two  «<U  ibrafcione  juai 
described  can  be  comp  a red  with  the  old  feutryl  calibration,  Columns 
l  and  2  of  Table  1ft  are  f rum  Column  1  of  Tables  U  and  VH  respectively, 
That  i a f  the  calibrated  pressure  la  given  at  1ft  yard  intervals  in 
Table  1ft,  Co  l  uron  3  represents  the  previous  calibration  for  tetryij 
an  interpolation  was  performed  on  the  previously  reported  results, 

The  diffnrenue  at  gero  gap  is  due  to  the  way  the  nominal  part  of 
the  calibration  was  dona,  In  both  the  new  calibrations*  an  ENV  gage 
was  used  to  measure  and  estimate  the  particle  velocity  at  the  explo» 
siv-  i'MMA  interface,  In  the  old  tetryl  oal  i  brat  ion,  the  interface 
value  was  computed,  using  CXJ  parameters  for  the  explosive  gases, 

This  resulted  in  a  lower  interface  pressure  than  that  which  is  orodueed 
in  the  explosive  ewnt,  If  we  had  used  values  for  the  parameters 
describing  the  state  in  the  explosive  gases  wluoh  is  appropriate 
to  the  von  Neumann  spike,  our  previous  resviUs  would  have  been  closer 
to  the  present  ones,  All  this  is  to  say  that  the  new  nominal  values 
are  better  approximations  than  those  published  previously. 

At  30  cards,  the  three  calibrations  are  remarkably  close/  the 
new  value  for  tetryl  is  slightly  smaller  than  the  old  value  and  the 
value  for  pentolite  is  slightly  greater,  as  expected,  At  the  end  of 
the  nominal  range  (3ft  cards),  the  new  value  for  to try. I  is  about 
1  kbar  larger  than  the  old,  The  difference  increases  to  a  maximum 
value  of  7  kbar  at  140  cards,  then  decreases  to  2,6  kbar  at  2oQ  cards. 
Thereafter  it  continues  to  decrease  to  0.4  kbar  at  320  earth*.  Prom 
330-400  cards,  the  two  values  are  the  same  *0.2  kbar. 

The  maximum  difference  in  pressure  (26  kbar  at  x  «  0)  amounts  to 
about  17*  change/  the  maximum  percentage  difference  is  20*  at  140 
caida.  The  uncertainty  in  the  earlier  calibration  was  estimated  s 
'  l  u  *  /  in  the  present  one  it.  is  *ib*.  The  maximum  change  is  there  tore 
about  the  sum  of  the  two  uncertainties, 

Respite  the  fact  that  the  new  results  for  tetryl  are  in  same 
ranges  significantly  larger  titan  the  early  ones,  the  two  (old  and 
now)  P  va  x  curves  do  not  cross.  Therefore  the  relative  ranking  of 
explosives  by  gap  sensitivity  will  not  be  changed  although  the  $0* 
pressures  (Pg)  will  be,  in  soma  oases.  The  same  situation  also  exists 
for  the  old  and  new  calibrations  with  pentolite. 

The  results  for  pentolite  are  consistently  (but  not  always 
significantly)  greater  than  those  of  tae  tetryl  calibration  except 
at  large  attenuations.  For  practical  purposes,  we  can  use  the  same 
calibration  (P  va  x)  curve  for  both  tetryl  and  pentolite  at  x  -s  200 
cards , 
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TabU  19 

COMPARISON  OF  PRKSSURR  CALIBRATION 
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7,1 

34  0 

6,7 

350 

6.3 

360 

6.0 

370 

5,7 

3  80 

5,5 

390 

5.2 

400 

5.0 

213,1 

13  4.8 

142,0 

130,3 

112,3 

m.o 

97,2 

96.3 

87,7 

83,0 

80,7 

76.3 

74,9 

69,3 

69,8 

63.4 

65,2 

38.2 

61.0 

33,7 

57,1 

49.5 

53,4 

45.6 

30,1 

41.9 

46.9 

3  8.1 

44.0 

34 ,4 

38.9 

30.8 

34.0 

27.4 

29,9 

24,3 

26 . 3 

21,6 

23.3 

19.3 

20,7 

17.3 

18.6 

15.7 

16.9 

14,3 

15.3 

13,1 

14.0 

12.1 

12.9 

11,2 

11.8 

10.3 

10.9 

9.6 

10. 1 

9,0 

9.3 

8.4 

8.7 

8.0 

8.1 

7.5 

7.6 

7.2 

7.1 

6.9 

6.7 

6.6 

6.3 

6.3 

6.0 

6.1 

5.7 

5.9 

5.4 

5.6 

5.2 

5.4 

5.0 
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APPENDIX  A 

ESTIMATION  OF  ERRORS  FOR  MEASUREMENT  OF  PARTICLE  VELOCITY 
WITH  THE  ELECTROMAGNETIC  GAGE 


The  basic  equation  for  the  electromagnetic  gage  is 


u  =  104  v/H l 


(Al) 


where  u  is  the  particle  velocity  (mm/Msec) ,  v  is  the  output  of  the 
gage  (volts),  H  is  the  magnetic  field  strength  (gauss)  and  i  is  the 
length  of  the  base  of  the  gage  (mm) .  The  quantity  l  is  measured 
directly  with  a  micrometer  so  that  an  estimate  of  error  is  straight¬ 
forward.  H  is  measured  with  a  gaussmeter  prior  to  the  shot  as  will 
be  described  later.  The  quantity  v  is  measured  indirectly  and  can 
be  in  error  because  of  several  sources.  It  is  measured  by  using  an 
oscilloscope  which  must  be  calibrated  before  each  shot.  A  reference 
line  is  placed  on  the  film  by  triggering  the  scope  manually,  see 
Figure  Al .  A  voltage,  E,  is  then  applied  to  the  input  to  the  scope 
and  it  is  manually  triggered  again,  giving  reference  line  2.  Then 
the  experiment  is  performed,  giving  a  trace  —  see  Figure  Al  for  an 
illustrative  example.  When  the  film  record  is  read,  the  separation 
of  the  two  reference  lines  is  measured.  Call  this  distance  x  as 
shown  in  the  figure.  Next,  the  trace  is  read;  this  amounts  to  mea¬ 
suring  the  vertical  distance  from  the  base  line  --  the  dashed  line 
in  Figure  Al  —  to  a  point  on  the  trace.  For  a  particular  point  on 
the  trace,  call  this  displacement  y.  We  assume  that  the  scope  is 
linear  so  that  displacement  of  the  trace  is  proportional  to  the 
voltage,  or 


V]_  =  E  y  /x. 


(A2) 


The  voltage  v^  is  not  that  generated  by  the  gage;  it  is  the  output  of 
the  circuit  shown  in  Figure  A2.  The  foil  is  represented  as  a  voltage 
generator  having  negligible  resistance  with  a  resistance  r  in  series. 
The  reason  for  incorporating  this  resistance  is  given  el ° "’"here . 

The  load  resistor,  R,  is  that  of  the  terminating  resistor  at  the 
oscilloscope.  From  the  diagram  it  is  evident  that 


v  =  vj_  (R  +  r)  /R  =  E  y  (R  +  r)  /  x  R  (A3) 
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REFERENCE  LINE  2 


REFERENCE  LINE  1 

FIG .  A1  IDEALIZED  OSCILLOSCOPE  RECORD 
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Thus  the  particle  velocity  is  given  by 

u  =  104  E  y  (R  +  r)/(H£  x  R)  . 


In  order  to  estimate  the  error  in  u,  Equation 
tiated  logarithmetically ,  giving 


Au  _  A£  4.  AE  4.  Ax  +  A/£  +  ( AR+Ar)  +  AR 

u  SL  E  x  y  (R+r)  R 


where  all  negative  signs  have  been  changed  to  positive.  We  must 
estimate  the  uncertainty  in  each  quantity,  Ae,  AJt,  etc.,  so  that 
their  cumulative  effect  on  the  particle  velocity  can  be  calculated. 
This  is  the  simplest  way  to  estimate  the  error.  A  more  sophisticated 
method  (which  gives  more  optimistic  results)  is  to  compute  Au/u  by 
taking  the  root  mean  square  of  the  values  of  the  terms  on  the  right 
hand  side  of  Equation  (A5).  Note  that  the  result  will  be  a  relative 
value,  a  nondimensional  number.  It  will  be  more  convenient  to  deal 
with  percentages  in  the  following  so  that  each  side  of  Equation  (A5) 
is  multiplied  by  100%. 

The  first  term  involves  l,  the  length  of  the  base  of  the  EMV  gage. 
For  most  of  our  work,  this  length  is  determined  by  the  thickness  of 
a  piece  of  PMMA  or  explosive  around  which  the  foil  is  wrapped  part 
way.  The  thickness  of  the  test  material  varies  from  2.5  to  10.0  mm 
and  is  measured  using  a  micrometer.  Assuming  these  measurements  are 
accurate  to  ±0.013  mm  (0.0005  inch),  the  accuracy  is  ±0.48%  for  the 
short  pickup  and  ±0.12%  for  the  longest  pickup.  The  reason  for  using 
short  pickups  is  to  minimize  the  effects  of  curvature  on  the  rise 
time  of  the  signal.  That  is,  when  we  place  a  gage  within  50  mm  of  the 
point  of  initiation  of  an  explosive,  we  use  a  short  pickup.  Some  of 
our  close-in  measurements  may  be  in  error  due  to  the  inaccurate  mea¬ 
surements  of  the  base  length  of  the  gage.  Note,  however,  under  these 
circumstances,  the  particle  velocity  is  large,  so  that  other  sources 
of  errors  are  at  their  minimum. 

The  next  term,  AE/E,  is  the  error  in  measuring  the  volt-  ge  applied 
to  the  scope  for  calibration.  This  voltage  is  supplied  by  a  small 
DC  power  supply.  It  is  measured  with  a  Model  7050  Fairchild  digital 
multimeter.  For  temperatures  between  15  and  35 °C,  the  accuracy  of  the 
multimeter  is  claimed  to  be  ±0.1%  of  the  reading,  ±1  digit.  We  use 
the  1.5  volt  scale  so  that  one  digit  is  10 “3  volts.  Thus  the  %  error 
is  ±0.1  (1.0  +  E)  /  E;  this  function  is  plotted  in  Figure  A3.  Our 
estimated  error  obviously  becomes  more  troublesome  when  we  have  to 
measure  small  voltages.  Our  lowest  voltage  is  about  0.15  volts  for 
which  we  can  expect  an  error  of  about  1%.  The  error  due  to  this 


(A4) 


(A4 )  is  differen- 
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source  is  the  same  for  all  points  read  off  any  one  record.  It  will 
differ,  of  course,  from  one  record  to  the  next  to  the  extent  to  which 
we  change  the  calibrating  voltage. 

The  quantity,  x,  is  a  measure  of  the  displacement  of  the  oscillo¬ 
scope  beam  due  to  the  calibrating  voltage,  E.  This  displacement  is 
read  with  a  Universal  Telereader.  It  is  doubtful  if  the  measuring 
instrument  is  limiting  our  accuracy  in  this  case.  Rather,  our  accu¬ 
racy  is  limited  by  the  ability  of  the  operator  to  place  a  crosswire 
accurately  on  the  image  of  the  trace.  These  traces,  which  are  placed 
on  the  film  prior  to  a  shot,  are  usually  narrow  and  distinct  and  are 
therefore  relatively  easy  to  measure.  The  operator  claims  that  a 
reading  can  be  repeated  to  within  ±10  counts.  The  distance  being  mea¬ 
sured  on  the  film  is  usually  about  20  mm,  or  1600  counts  on  the 
Telereader.  If  both  lines  are  read  to  the  same  precision,  this 
distance  is  uncertain  to  ±20  counts,  or  1.25%.  This  error  is  the  same 
for  any  one  record.  It  does  not  change  much  from  record  to  record  even 
though  the  particle  velocity  may  vary  by  a  factor  of  20.  This  is  due 
to  the  fact  that  we  change  the  gain  of  the  amplifiers  of  the  scope  to 
give  about  20  mm  displacement  each  time  we  calibrate  the  scope.  We 
assume,  of  course,  that  the  amplifiers  are  linear. 

The  next  term  involves  y,  the  quantity  which  represents  the 
displacement  of  a  point  on  the  dynamic  record  from  the  base  line,  see 
Figure  Al.  In  the  usual  case,  we  measure  a  displacement  of  about 
10  mm  on  the  film,  or  about  800  count,  on  the  Telereader  scale.  Again 
we  assume  that  the  Telereader  is  so  accurate  that  it  does  not  limit 
our  accuracy.  The  crosswires  can  be  reset  on  a  trace,  and  on  the 
baseline,  to  within  about  ±10  counts.  Thus  the  displacement  of  the 
trace  is  uncertain  to  ±20  counts,  or  about  2.6%.  As  in  the  case  of 
the  measurement  of  x  as  discussed  above,  the  error  due  to  the  measure¬ 
ment  of  y  does  not  change  from  record  to  record.  The  reason  is  the 
same,  we  change  the  amplification  of  the  scope  so  that  the  trace  is 
displaced  about  the  same  in  each  shot.  Because  we  are  usually  mea¬ 
suring  a  decaying  signal,  the  error  does  change  some  over  each  indi¬ 
vidual  trace. 

Reported  results  never  depend  on  the  reading  of  a  single  point 
on  the  trace  of  a  iccord.  Many  points  are  read  so  that  a  plot  of  u 
vs  time  can  be  constructed.  A  subset  of  the  points  may  be  fitted  to 
a  straight  line  so  that  the  value  of  u  can  be  obtained  at  zero  time 
(an  extrapolation  back  through  the  rise  time  of  the  signal) .  In  other 
cases,  lines  are  drawn  through  the  data  to  locate  more  or  less  distinct 
changes  of  slope  of  the  trace.  This  graphical  treatment  has  been  used 
to  determine  the  Chapman-Jouguet  particle  velocity  and  the  reaction 
time  in  explosives.  Analysis  of  the  errors  under  these  more  compli¬ 
cated  circumstances  has  not  been  made.  The  error  reported  here 
applies  to  the  reading  of  only  a  single  point  on  the  trace  of  an 
oscilloscope  record. 

We  use  a  Rawson-Lush  Model  824  rotating  coil  gaussmeter  to  measure 
the  field,  H.  The  accuracy  of  this  instrument  is  0.1%  or  0.5  gauss, 
whichever  is  larger.  We  work  with  fields  whose  strength  vary  from 
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400  to  1000  gausses,  so  that  our  accuracy  varies  from  0.12%  to  0.1%. 
Errors  of  this  magnitude  are  negligible.  Our  method  of  operation  may 
introduce  other  errors  in  the  field,  however.  The  magnet  is  energized 
and  the  field  is  measured  before  the  experiment  is  completely  assem¬ 
bled.  The  current  to  the  magnet  is  measured  on  a  0-50  ampere  ammeter 
calibrated  to  0.25%  of  full  scale.  The  current  is  then  reduced  to 
zero,  the  leads  to  the  magnet  are  disconnected  for  reasons  of  safety, 
and  the  experiment  including  uhe  explosive,  is  put  into  place.  Imme¬ 
diately  before  firing,  the  current  to  the  magnet  is  restored  to  the 
same  value  as  before.  Firing  takes  place  within  3  to  5  minutes  of 
the  time  that  the  current  is  restored.  Drift  in  the  output  of  the 
power  supply  during  this  short  time  is  negligible.*  If  the  magnet 
assembly  was  not  disturbed  during  the  final  assembly  of  the  experiment 
and  if  the  ammeter  readings  were  the  same,  we  should  have  the  same 
field  as  that  during  the  measurement  with  the  gaussmeter.  Reading 
the  ammeter  may  introduce  an  error  in  the  field  strength.  The  meter 
has  '  TM.*'-:or  scale  which  helps  eliminate  parallax  during  reading.  On 
repeated  trials  to  check  on  this  possibility,  the  field  is  returned 
to  within  1.0  gauss  of  the  original  reading.  This  gives  an  error  of 
0.25%  at  400  gauss  and  0.1%  at  1000  gauss.  We  hope  that  in  routine 
operation  that  we  do  this  well. 

Both  of  the  resistors  in  the  circuit,  see  Figure  A1 ,  have  nominal 
values  of  50  ohms.  They  are  used  to  terminate  the  signal  cable  so 
that  the  reflection  of  electrical  signals  will  be  minimized.  The 
precision  of  the  resistor,  R,  which  is  located  at  the  oscilloscope, 
is  ±0.5%.  Hence  the  value  of  AR/R  is  also  ±0.5%.  The  resistor,  r, 
is  an  inexpensive  carbon  resistor  which  is  destroyed  in  the  experiment 
Its  resistance  is  measured  on  a  bridge  to  a  precision  of  ±0.25%  prior 
to  assemblying  the  shot.  Converting  the  percentages  given  above  to 
increments  of  resistances  gives  0.38%  for  the  term  (AR  +  Ar)/(R  +  r) . 
Hence  we  have  a  total  error  from  the  resistors  of  0.88%.  This  value 
can  be  reduced  by  a  careful  measurement  of  the  resistance  at  the 
oscilloscope  (resistor  R) .  Note  that  in  the  above  that  the  resistance 
of  the  cable  has  been  neglected. 

The  errors  are  summarized  below  for  three  values  of  u.  This 
illustrates  how  the  error  in  the  three  digit  voltmeter  changes  the 
error  as  u  is  changed.  The  error  in  y  is  included  in  the  table  so 
that  the  results  do  not  apply,  for  example,  to  the  determination  of 
a  peak  value  of  u  which  usually  depends  on  fitting  many  values  of  u 
to  a  function  for  extrapolation  to  zero  time. 

In  the  analysis  of  error  given  above,  we  assumed  that  experimental 
arrangements  were  reproducible.  There  may  be  variations  in  the 


*The  power  supply  is  an  Electronic  Measurements  'del  SCR  120-20, 
three  phase,  constant  current  instrument. 
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SUMMARY  OF  ERRORS 


u,  mm/psec 

2.0 

1.0 

0.1 

M/l 

0.50 

0.5 

0.5 

AE/E 

0.14 

0.24 

0.75 

A  x/x 

1.26 

1.26 

1.26 

Ay/y 

2.6 

2.6 

2.6 

(AR+Ar)/  (R+r) 

0.38 

0.38 

0.38 

AR/R 

0.5 

0.5 

0.5 

AH/H 

0.25 

0.25 

0.25 

Sum  of  Errors 

5.63 

5.73 

6.24 

density  of  the  PMMA  from  sample  to  sample.  The  same  is  true  of  the 
pentolite  pellets.  The  PMMA  samples  are  machined  in  a  lathe  so  that 
they  are  at  least  in  the  correct  geometrical  shape.  The  pressed 
pentolite  pellets  are  somewhat  barrel  shaped  —  the  diameters  of  the 
midplane  of  the  pellets  are  about  0.003  inch  greater  than  the  diam¬ 
eter  of  their  faces.  The  faces  are  bulged  so  that  on  the  axis  the 
pellets  are  0.005  to  0.006  inch  thicker  than  they  are  at  their  edges. 
Thus  when  two  pellets  are  stacked,  it  is  possible  to  have  a  gap 
0.012  inch  thick  over  part  of  the  area  between  them.  This  interferes 
with  the  transfer  of  detonation  from  one  pellet  to  another,  introduces 
a  tilt  in  the  detonation  wave,  and  distorts  the  Taylor  wave.  At  the 
explosive/PMMA  interface,  the  gap  may  be  0.006  inch  over  part  of  the 
surface  (assuming  the  PMMA  is  flat  on  the  end) .  Again  this  introduces 
tilt  in  the  shock  wave  in  the  PMMA  and  distorts  the  rarefaction  wave 
which  follows  the  shock  front. 

The  shape  of  tetryl  pellets  are  closer  to  that  of  a  cylinder. 

They  may  be  as  much  as  0.002  inch  thicker  on  their  axes  than  at  their 
edges.  Their  diameters  may  be  0.002  inch  too  great  at  their  midplanes. 

Our  calibration  work  might  have  been  easier  if  we  had  machined  the 
pressed  pellets  into  true  cylindrical  forms.  Then,  however,  the 
calibration  would  not  have  corresponded  to  the  actual  situation  in 
the  gap  test. 
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